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Abstract. In recent years, quantum phase transitions have attracted the interest of 
both theorists and experimentalists in condensed matter physics. These transitions, 
which are accessed at zero temperature by variation of a non-thermal control 
parameter, can influence the behavior of electronic systems over a wide range of the 
phase diagram. Quantum phase transitions occur as a result of competing ground state 
phases. The cuprate superconductors which can be tuned from a Mott insulating to 
a d-wave superconducting phase by carrier doping are a paradigmatic example. This 
review introduces important concepts of phase transitions and discusses the interplay 
of quantum and classical fluctuations near criticality. The main part of the article 
is devoted to bulk quantum phase transitions in condensed matter systems. Several 
classes of transitions will be briefly reviewed, pointing out, e.g., conceptual differences 
between ordering transitions in metallic and insulating systems. An interesting 
separate class of transitions are boundary phase transitions where only degrees of 
freedom of a subsystem become critical; this will be illustrated in a few examples. The 
article is aimed on bridging the gap between high-level theoretical presentations and 
research papers specialized in certain classes of materials. It will give an overview on 
a variety of different quantum transitions, critically discuss open theoretical questions, 
and frequently make contact with recent experiments in condensed matter physics. 
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1. Introduction 

Phase transitions play an essential role in nature. Everyday examples include the boiling 
of water or the melting of ice, more complicated is the transition of a metal into the 
superconducting state upon lowering the temperature. The universe itself is thought to 
have passed through several phase transitions as the high-temperature plasma formed 
by the Big Bang cooled to form the world as we know it today. 

Phase transitions occur upon variation of an external control parameter; their 
common characteristics is a qualitative change in the system properties. The phase 
transitions mentioned so far occur at finite temperature; here macroscopic order (e.g. 
the crystal structure in the case of melting) is destroyed by thermal fluctuations. 
During recent years, a different class of phase transitions has attracted the attention of 
physicists, namely transitions taking place at zero temperature. A non-thermal control 
parameter such as pressure, magnetic field, or chemical composition, is varied to access 
the transition point. There, order is destroyed solely by quantum fluctuations which are 
rooted in the Heisenberg uncertainty principle. 

Quantum phase transitions pi Ej have become a topic of vivid interest in 
current condensed matter physics. At first glance it might appear that the study of 
such special points in the phase diagram is a marginal problem of interest only to 
specialists, as such transitions occur at only one special value of a control parameter 
at the experimentally impossible temperature of absolute zero. However, experimental 
and theoretical developments in the last decades have clearly established the contrary. 
They have made clear that the presence of such zero-temperature quantum critical 
points holds the key to so-far unsolved puzzles in many condensed matter systems. 
Examples include rare-earth magnetic insulators heavy-fermion compounds [SJ Ej 
high-temperature superconductors [HJ E] , and two-dimensional electron gases [TJ [TU] . 

As we will see below, quantum critical behavior, arising from the peculiar excitation 
spectrum of the quantum critical ground state, can influence measurable quantities over 
a wide range of the phase diagram. The physical properties of the quantum fluctuations, 
which can destroy long-range order at absolute zero, are quite distinct from those of the 
thermal fluctuations responsible for traditional, finite-temperature phase transitions. In 
particular, the quantum system is described by a complex-valued wavefunction, and the 
dynamics of its phase near the quantum critical point requires novel theories which have 
no analogue in the traditional framework of phase transitions. 

This article is intended as primarily non-technical introduction into the field of 
quantum phase transitions, and could serve as a reference for both theorists and 
experimentalists interested in the field. In Sec. El we start with summarizing the 
basics of finite-temperature phase transitions, extend those general concepts to T = 
highlighting the interplay between classical and quantum fluctuations near a quantum 
critical point, and illustrate the correspondence between quantum transitions in d and 
classical transitions in D = d + z dimensions, where z is the dynamic critical exponent. 
In addition, we briefly discuss some aspects of the theoretical description of phase 
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transitions, like order parameter field theories and their renormalization group analysis. 
Sec. El is devoted to the concrete description of several classes of bulk quantum phase 
transitions, where "bulk" refers to the fact that the whole system becomes critical. The 
discussion includes conventional transitions involving spin or charge order as well as 
metal-insulator and superconductor-insulator transitions; it is intended to summarize 
some of the settled and open issues in the field, and to make contact with experiments. 
For more technical details the reader is referred to the references. In Sec. Ejwe introduce 
so-called boundary quantum phase transitions - those are transitions where only degrees 
of freedom of a subsystem become critical. This concept is illustrated using examples 
of quantum impurity problems, which are of current interest both in correlated bulk 
materials and in mesoscopic physics. We end in Sec. 03 with a brief summary and 
outlook. To shorten notations, we will often employ units such that % = ks = 1. 

2. Quantum phase transitions 

This section will give a general introduction into both classical and quantum phase 
transitions, and point out similarities and important conceptual differences between the 
two. 

2.1. Basic concepts of phase transitions 

We start out with briefly collecting the basic concepts of phase transitions and critical 
behavior fTJ E] which are necessary for the later discussions. Phase transitions 
are traditionally classified into first-order and continuous transitions. At first-order 
transitions the two phases co-exist at the transition temperature - examples are ice and 
water at 0° C, or water and steam at 100° C. In contrast, at continuous transitions the 
two phases do not co-exist. An important example is the ferromagnetic transition of iron 
at 770° C, above which the magnetic moment vanishes. This phase transition occurs at 
a point where thermal fluctuations destroy the regular ordering of magnetic moments 
- this happens continuously in the sense that the magnetization vanishes continuously 
when approaching the transition from below. The transition point of a continuous phase 
transition is also called critical point. The study of phase transitions, continuous phase 
transitions in particular, has been one of the most fertile branches of theoretical physics 
in the last decades. 

In the following we concentrate on systems near a continuous phase transition. 
Such a transition can usually be characterized by an order parameter - this is a 
thermodynamic quantity that is zero in one phase (the disordered) and non-zero and 
non-unique in the other (the ordered) phase. Very often the choice of an order parameter 
for a particular transition is obvious as, e.g., for the ferromagnetic transition where 
the total magnetization is an order parameter. However, in some cases finding an 
appropriate order parameter is complicated and still a matter of debate, e.g., for the 
interaction-driven metal-insulator transition in electronic systems (the Mott transition 
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While the thermodynamic average of the order parameter is zero in the disordered 
phase, its fluctuations are non-zero. If the critical point is approached, the spatial 
correlations of the order parameter fluctuations become long-ranged. Close to the critical 
point their typical length scale, the correlation length £, diverges as 

e«i*r (i) 

where v is the correlation length critical exponent and t is some dimensionless measure 
of the distance from the critical point. If the transition occurs at a non-zero temperature 
T c , it can be defined as t = \T — T c \/T c . In addition to the long-range correlations in 
space there are analogous long-range correlations of the order parameter fluctuations 
in time. The typical time scale for a decay of the fluctuations is the correlation (or 
equilibration) time r c . As the critical point is approached the correlation time diverges 
as 

r c oc i z oc \t\- vz (2) 

where z is the dynamic critical exponent. Close to the critical point there is no 
characteristic length scale other than £ and no characteristic time scale other than 
t c . (Note that a microscopic cutoff scale must be present to explain non-trivial critical 
behavior, for details see, e.g., Goldenfeld [Ej. In a solid such a scale is, e.g., the lattice 
spacing.) 

The divergencies (0) and (J2J) are responsible for the so-called critical phenomena. 
At the phase transition point, correlation length and time are infinite, fluctuations 
occur on all length and time scales, and the system is said to be scale-invariant. As a 
consequence, all observables depend via power laws on the external parameters. The 
set of corresponding exponents - called critical exponents - completely characterizes the 
critical behavior near a particular phase transition. 

Let us illustrate the important concept of scaling in more detail. To be specific, 
consider a classical ferromagnet with the order parameter being the magnetization M(r). 
External parameters are the reduced temperature t = \T — T c \/T c and the external 
magnetic field B conjugate to the order parameter. Close to the critical point the 
correlation length is the only relevant length scale, therefore the physical properties 
must be unchanged if we rescale all lengths in the system by a common factor, and at 
the same time adjust the external parameters in such a way that the correlation length 
retains its old value. This gives rise to the homogeneity relation for the singular part of 
the free energy density, 

f(t,B) = b- d f(tbV»,Bby»). (3) 

Here ys is another critical exponent. The scale factor b is an arbitrary positive number. 
Analogous homogeneity relations for other thermodynamic quantities can be obtained 
by differentiating /. The homogeneity law (JHJ) was first obtained phenomenologically by 
Widom within the framework of the renormalization group theory it can be 
derived from first principles. 
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Table 1. Commonly used critical exponents for magnets, where the order parameter 
is the magnetization m, and the conjugate field is a magnetic field B. t denotes the 
distance from the critical point and d is the space dimensionality. (The exponent ub 
defined in J3J) is related to 8 by ys = dS/(l + 5).) 
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In addition to the critical exponents u, db and z defined above, a number of other 
exponents is in common use. They describe the dependence of the order parameter 
and its correlations on the distance from the critical point and on the field conjugate 
to the order parameter. The definitions of the most commonly used critical exponents 
are summarized in Table H Note that not all the exponents defined in Table are 
independent from each other. The four thermodynamic exponents a, (3, 7, S can all be 
obtained from the free energy Q which contains only two independent exponents. They 
are therefore connected by the so-called scaling relations 

2- a = 2/3 + 7, 2-a = p{8+\) . (4) 

Analogously, the exponents of the correlation length and correlation function are 
connected by two so-called hyperscaling relations 

2-a = dv, 7 = (2 -77)1/. (5) 

(Hyperscaling relations are violated in theories with mean-field critical behavior due 
to the presence of a dangerously irrelevant variable, see Sec. 12.61 ) Since statics and 
dynamics decouple in classical statistics (Sec. 12.3(1 the dynamic exponent z is completely 
independent from all the others. 

One of the most remarkable features of continuous phase transitions is universality, 
i.e., the fact that the critical exponents are the same for entire classes of phase 
transitions which may occur in very different physical systems. These universality 
classes are determined only by the symmetries of the order parameter and by the 
space dimensionality of the system. This implies that the critical exponents of a 
phase transition occurring in nature can be determined exactly (at least in principle) 
by investigating any simple model system belonging to the same universality class. The 
mechanism behind universality is again the divergence of the correlation length. Close 
to the critical point the system effectively averages over large volumes rendering the 
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microscopic details of the Hamiltonian unimportant. 

2.2. Quantum mechanics and the vicinity of the critical point 

The question of to what extent quantum mechanics is important for understanding 
a continuous phase transition has at least two aspects. On the one hand, quantum 
mechanics can be essential to understand the existence of the ordered phase, (e.g., 
superconductivity) - this depends on the particular transition considered. On the 
other hand, one may ask whether quantum mechanics influences the asymptotic critical 
behavior. For this discussion we have to compare two energy scales, namely hu> c , which 
is the typical energy of long-distance order parameter fluctuations, and the thermal 
energy k B T . We have seen in the preceeding section that the typical time scale r c of 
the fluctuations diverges as a continuous transition is approached. Correspondingly, the 
typical frequency scale uj c goes to zero and with it the typical energy scale 

hu c oc \t\ uz . (6) 

Quantum mechanics will be important as long as this typical energy scale is larger than 
the thermal energy k B T; on the other hand, for %uj c <C k B T a purely classical description 
can be applied to the order parameter fluctuations. In other words, the character of 
the order parameter fluctuations crosses over from quantum to classical when hu c falls 
below k B T. 

Now, for any transition occurring at some finite temperature T c quantum mechanics 
will become unimportant for \t\ < T^ uz , in other words, the critical behavior 
asymptotically close to the transition is entirely classical. This justifies to call all finite- 
temperature phase transitions "classical". Quantum mechanics can still be important 
on microscopic scales, but classical thermal fluctuations dominate on the macroscopic 
scales that control the critical behavior. If, however, the transition occurs at zero 
temperature as a function of a non-thermal parameter r like pressure or magnetic field, 
the behavior is always dominated by quantum fluctuations. Consequently, transitions 
at zero temperature are called "quantum" phase transitions. 

The interplay of classical and quantum fluctuations leads to an interesting phase 
diagram in the vicinity of the quantum critical point. Two cases need to be distinguished, 
depending on whether long-range order can exist at finite temperatures. 

Fig. UK describes the situation where order only exists at T = 0, this is the case, e.g., 
in two-dimensional magnets with SU(2) symmetry where order at finite T is forbidden 
by the Mermin- Wagner theorem. In this case there will be no true phase transition in 
any real experiment carried out at finite temperature. However, the finite-T behavior 
is characterized by three very different regimes, separated by crossovers, depending on 
whether the behavior is dominated by thermal or quantum fluctuations of the order 
parameter. In the thermally disordered region the long-range order is destroyed mainly 
by thermal order parameter fluctuations. In contrast, in the quantum disordered region 
the physics is dominated by quantum fluctuations, the system essentially looks as in 
its quantum disordered ground state at r > r c . In between is the so-called quantum 
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Figure 1. Schematic phase diagrams in the vicinity of a quantum critical point 
(QCP). The horizontal axis represents the control parameter r used to tunc the system 
through the quantum phase transition, the vertical axis is the temperature T. a) 
Order is only present at zero temperature. The dashed lines indicate the boundaries of 
the quantum critical region where the leading critical singularities can be observed; 
these crossover lines are given by fc^T oc \r — r c \ vz . b) Order can also exist at 
finite temperature. The solid line marks the finite-temperature boundary between 
the ordered and disordered phases. Close to this line, the critical behavior is classical. 



critical region [T^j, where both types of fluctuations are important. It is located near the 
critical parameter value r = r c at comparatively high (!) temperatures. Its boundaries 
are determined by the condition ksT > hu c oc \r — r c \ vz : the system "looks critical" 
with respect to the tuning parameter r, but is driven away from criticality by thermal 
fluctuations. Thus, the physics in the quantum critical region is controlled by the 
thermal excitations of the quantum critical ground state, whose main characteristics is 
the absence of conventional quasiparticle-like excitations. This causes unusual finite- 
temperature properties in the quantum critical region, such as unconventional power 
laws, non- Fermi liquid behavior etc. Universal behavior is only observable in the vicinity 
of the quantum critical point, i.e., when the correlation length is much larger than 
microscopic length scales. Quantum critical behavior is thus cut off at high temperatures 
when ksT exceeds characteristic microscopic energy scales of the problem - in magnets 
this cutoff is, e.g., set by the typical exchange energy. 

If order also exists at finite temperatures, Fig. the phase diagram is even richer. 
Here, a real phase transition is encountered upon variation of r at low T; the quantum 
critical point can be viewed as the endpoint of a line of finite-temperature transitions. 
As discussed above, classical fluctuations will dominate in the vicinity of the finite-T 
phase boundary, but this region becomes narrower with decreasing temperature, such 
that it might even be unobservable in a low-T experiment. The fascinating quantum 
critical region is again at finite temperatures above the quantum critical point. 

A quantum critical point can be generically approached in two different ways: 
r — > r c at T = or T — > at r = r c . The power-law behavior of physical observables in 
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both cases can often be related. Let us discuss this idea by looking at the entropy S. It 
goes to zero at the quantum critical point (exceptions are impurity transitions discussed 
in Sec. EJ), but its derivatives are singular. The specific heat C = TdS/dT is expected 
to show power-law behavior, as does the quantity B = dS/dr. Using scaling arguments 
(see below) one can now analyze the ratio of the singular parts of B and C: the scaling 
dimensions of T and S cancel, and therefore B/C scales as the inverse of the tuning 
parameter r. Thus, one obtains a universal divergence in the low-temperature limit, 
B/C oc \r - r c \~ l ] similarly, B/C oc T~ x l vz at r = r c [17 a . Note that B/C does not 
diverge at a finite-temperature phase transition. At a pressure-tuned phase transition 
r = p, then B measures the thermal expansion, and B/C is the so-called Griineisen 
parameter. As will be discussed in Sec. 12.61 the scaling argument presented here can be 
invalid above the upper-critical dimension. 

2.3. Quantum- classical mapping and scaling 

To gain a deeper understanding of the relation between classical and quantum behavior, 
and the possible quantum-classical crossover, we have to recall general features from 
quantum statistical mechanics. 

The starting point for the derivation of thermodynamic properties is the partition 
function 

Z = Tre -H/k B T (?) 

where H = ifkm + -f^pot is the Hamiltonian characterizing the system. In a classical 
system, the kinetic and potential part of H commute, thus Z factorizes, Z = Z kin Z pot , 
indicating that in a classical system statics and dynamics decouple. The kinetic 
contribution to the free energy will usually not display any singularities, since it derives 
from the product of simple Gaussian integrals. Therefore one can study classical 
phase transitions using effective time-independent theories, which naturally live in d 
dimensions. 

In contrast, in a quantum problem the kinetic and potential parts of H in general 
do not commute, the quantum mechanical partition function does not factorize, which 
implies that statics and dynamics are always coupled. An order parameter field theory 
needs to be formulated in terms of space and time dependent fields. The canonical 
density operator e~ H ' kBT looks exactly like a time evolution operator in imaginary time r 
if one identifies 1/ksT = t = —iQ/h where G denotes the real time. Therefore it proves 
convenient to introduce an imaginary time direction into the system - formally this is 
done in the path integral representation of the partition function. At zero temperature 
the imaginary time acts similarly to an additional space dimension since the extension 
of the system in this direction is infinite. According to (J2J), time scales like the z-th 
power of a length. (Note that z—1 for many transitions in clean insulators, however, in 
general other values of z including fractional ones can occur.) The classical homogeneity 
law (jHJ) for the free energy density can now easily be adopted to the case of a quantum 
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phase transition. At zero temperature it reads 



f(t,B) = b- {d+z) f(tb 1/v ,Bb yB ) 



(8) 



where now t = \r — r c \/r c . This shows that a quantum phase transition in d space 
dimensions is related to a classical transition in (d+z) space dimensions. 

For a quantum system, it is now interesting to discuss the quantum to classical 
crossover upon approaching a finite-temperature phase transition - this turns out to 
be equivalent to a dimensional crossover. Under the quantum-classical mapping the 
temperature of the quantum problem maps onto the inverse length of the imaginary 
time axis. (Note that the temperature of the classical problem corresponds to a coupling 
constant in the quantum model within this mapping.) Close to the transition (Fig. ^d) 
the behavior is determined by the relation between the characteristic correlation time, 
r c , and the extension in imaginary time direction, (3 = 1/ksT. The crossover from 
quantum to classical behavior will occur when the correlation time r c becomes larger 
than (3 which is equivalent to the condition \t\ uz < ksT; in other words, once r c > (3 
the system realizes that it is effectively only <i-dimensional and not (cf+z)-dimensional. 
The corresponding crossover scaling is equivalent to finite size scaling in imaginary time 
direction. 

In contrast, when approaching a zero -temperature transition by lowering the 
temperature at r = r c , both 1/ksT and r c diverge, such that quantum effects are 
always important, and we get a truly (d + z)-dimensional system. Recognizing that 
quantum critical singularities can be cut off both by tuning r away from r c at T = 
and by raising T at r = r c , it is useful to generalize the homogeneity law (JHJ) to finite 
temperatures, 



The homogeneity law for the free energy immediately leads to scaling behavior of both 
static and dynamic observables. Consider an observable 0(k,u), e.g., a magnetic 
susceptibility, measured at wavevector k and frequency u>, then the existence of a single 
length scale £ and a single time scale u~ x = £ z implies 



where 0\ and O2 are different forms for the scaling function associated with the 
observable O, and do is the so-called scaling dimension of O. For simplicity we have 
set the external field B = and assumed that k is measured relative to the ordering 
wavevector. 

Precisely at the critical point, where the correlation length is infinite, the only length 
is set by the measurement wavevector k, analogously, the only energy is u, leading to 




(9) 



0(t,k,u,T) = £ d ° O^ouCTe) 




(10) 




(11) 



Similarly, for k = 0, but a finite temperatures, we obtain 




(12) 
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this is known as U u/T" scaling. It is important to note that the relations (jlU)) - (jl2)) . 
sometimes also called "naive scaling" , are only expected to be valid if the critical point 
satisfies hyperscaling properties, which is true below the upper-critical dimension, see 
Sec. EH 

The power of the quantum-classical analogy can be nicely demonstrated by 
considering a correlation function at the critical point. In a classical <i-dimensional 
system, correlations typically fall off with a power law in real space implying a 
momentum space behavior as 

G(k) oc k- 2+71d . (13) 

If we are now interested in a quantum phase transition with z = 1 which maps onto 
a (d + l)-dimensional classical theory with known properties, we can conclude that at 
T = 

G(k,iuj n ) oc [k 2 + oo 2 J- 2+ ^ 2 (14) 

where the Matsubara frequency uo n simply takes the role of an additional wavevector 
component. For real frequencies, the above translates into a retarded Green's function 
of the form 

G R (k, u) oc [k 2 + ^) 2 ](- 2 +^+i)/2 . (15) 

This examples illustrates the character of the excitation spectrum at a quantum critical 
point: G R does not show a conventional quasiparticle pole, but instead a branch cut for 
uj > k, corresponding to a critical continuum of excitations. This implies that modes 
become overdamped, and the system shows quantum relaxational dynamics 

Now, if a quantum transition seems to map generically onto a classical transition, 
what is different about quantum transitions? For a number of reasons it turns out that 
not all needed properties of a given quantum system can be obtained from a classical 
theory, the reasons being: 

(A) Theories for quantum systems can have ingredients which make them 
qualitatively different from their classical counterparts; examples are topological Berry 
phase terms and long-ranged effective interactions arising from soft modes. 

(B) Even if a quantum-classical mapping is possible, calculating real time dynamics 
of the quantum system often requires careful considerations and novel theories, as an 
approximation done on the imaginary time axis is in general not appropriate for real 
times 0. 

(C) Quenched disorder leads to an extreme anisotropy in space-time in the classical 
theory JTH] , see Sec. 13.31 

(D) The dynamic properties of the quantum system near a quantum critical 
point are characterized by a fundamental new time scale, the phase coherence time 
T0, which has no analogue at the corresponding classical transition [3]. Notably, 
diverges as T — > for all parameter values, i.e., the quantum system has perfect phase 
coherence even in the disordered phase - this seems to be peculiar considering that 
all correlations decay exponentially in the high-temperature phase of the corresponding 
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(c/+2;)-dimensional classical model. Technically, this peculiarity is related to the nature 
of the analytic continuation between imaginary and real times. In many models for 
quantum phase transitions, it is found that oc 1/T as T — > at the critical coupling, 
but diverges faster in the stable phases. 

Parenthetically, we note that "conventional" quantum criticality leads to a power 
law behavior in dynamic observables as function of energy and momentum, where (at 
T = 0) the singularity is cut off both by finite frequency and finite momentum (measured 
from the ordering wavevector). There are, however, interesting cases where the order 
parameter is not associated with a particular momentum, and then so-called "local" 
critical behavior may occur. Candidate transitions are those involving spin glass or 
topological order; interesting recent proposals involve the breakdown of Kondo screening 
in heavy fermion compounds (Sec. I3.2."3|) . 

2. 4- Ground state properties near a quantum phase transition 

It is useful to think about a quantum transition from the point of view of many-body 
eigenstates of the system. A quantum phase transition is a non-analyticity of the ground 
state properties of the system as function of the control parameter r. If this singularity 
arises from a simple level crossing in the many-body ground state, then we have a first- 
order quantum phase transition, without diverging correlations and associated critical 
singularities. A first-order quantum transition can also occur in a finite-size system. 
The situation is different for continuous transitions, where a higher-order singularity in 
the ground state energy occurs: here an infinite number of many-body eigenstates are 
involved, and the thermodynamic limit is required to obtain a sharp transition. For 
any finite-size system a continuous transition will be rounded into a crossover - this is 
nothing but an avoided level crossing in the ground state, which then becomes infinitely 
sharp in the thermodynamic limit. 

For numerical simulations of finite-size systems the choice of boundary conditions 
plays an important role: if the infinite system shows a first-order transition where the 
two ground states can be distinguished by a certain quantum number, then such a 
transition will also be rounded into a crossover for finite system size if the boundary 
conditions spoil the existence of this quantum number. This applies, e.g., to simulations 
using open boundary conditions which invalidate momentum or a certain parity as good 
quantum numbers. 

2. 5. Order parameters and order parameter field theories 

Near a phase transition the relevant physics is described by long- wavelength order 
parameter fluctuations. Therefore it appears appropriate to disregard microscopic 
processes and to work with a theory containing the order parameter fluctuations only. 
However, in itinerant electron systems it can be neccessary to include low-energy 
fermionic excitations explicitly in the critical theory jTH] because of the strong coupling 
between order parameter and fermionic particle-hole modes, see Sees. 13.2.31 and 13.2.41 
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Formally, an effective theory can be obtained in the path integral formulation from 
a microscopic model of interacting particles (e.g. fermions). To this end, the interaction 
terms are suitably decoupled via a Hubbard- Stratonovich transformation, and then the 
microscopic degrees of freedom can be integrated out, resulting in a theory of the order 
parameter fluctuations alone. Near a phase transition, the correlation length is large, 
and therefore the effective theory can be formulated in the continuum limit, with a 
proper ultraviolet cutoff set by the lattice spacing. In most cases the resulting theory 
will contain simple powers of the order parameter field and its gradients, and the form 
of such a theory can actually be derived based on symmetry arguments only. The 
form of the low-energy effective action will thus only depend on the dimensionality 
and symmetry of the underlying system and the symmetry of the order parameter - 
reflecting nothing but the celebrated universality of critical behavior. 

In the following we mention a few examples of models and corresponding quantum 
field theories; concrete applications will be given in Sec. El We will not discuss here 
the special theories appropriate for one-dimensional systems, the reader is referred to 
Ref. |3j for an introduction. 

2.5.1. Ising and rotor models. Quantum rotor models, with the Hamiltonian 

f^^E^-JE^i, (16) 

i (ij) 

are some of the simplest models showing quantum phase transitions. Here, the N- 
component operators with N > 2 represent the orientation of rotors on sites % of a 
regular rf-dimensional lattice, with nf = 1, and the L, are the corresponding angular 
momenta. The nearest-neighbor coupling in the second term prefers a ferromagnetic 
ordering, whereas the kinetic term tends to delocalize the individual rotors. Thus, 
j < 1 leads to an ordered state, whereas g > 1 places the system into a quantum- 
disordered phase [3]. A related model with N — 1 components is the quantum Ising 
model discussed in more detail in Sec. 12.81 

The order parameter describing the transition in the model (j!6|) is the ferromagnetic 
magnetization; a continuum theory can be derived from the lattice model in the path 
integral formulation by spatial coarse- graining, i.e., the variables fij are averaged over 
microscopic length scales. The magnitude of the continuum order parameter field, 0, can 
thus vary over a wide range, although the "magnitude" of the underlying microscopic 
object (a single rotor) is fixed. Expanding the resulting action into lowest-order field 
gradients, one arrives at the so-called 4 theory. This quantum field theory, described 
by the action 



Z= J V(J) a (x,T) exp(-<S) , 



S = I d d x I dr 



h/k B T 







(17) 



turns out to be a generic theory, which we will refer to frequently in this article, being 
appropriate for systems where the ^-component order parameter field, 
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non-conserved density, r is the bare "mass" of the order parameter fluctuations, used to 
tune the system across the quantum phase transition, and c is a velocity, go describes 
the important self-interactions of the <fi fluctuations, and we will comment on them in 
the next section. Note that the dynamic exponent is z — 1, i.e., space and time enter 
symmetrically. Although the order parameter fluctuations are bosonic in nature, the 
second-order time derivative in (j!7|) does not represent the dynamics of a canonical 
boson, but rather arises from a gradient expansion of the order parameter action. 

The model Eq. ()17j) obeys a O(N) symmetry, corresponding to rotations in order 
parameter space. N = 1 corresponds to a scalar or Ising order parameter; for N > 2 
the order parameter represents a vector for which one can distinguish longitudinal 
(amplitude) and transverse (phase or direction) fluctuations. For N = 3, the theory 
can describe spin degrees of freedom in an insulator; in a number of cases, however, 
Berry phases are important which we will explicitely discuss below. In writing down 
(JT7J) we have assumed to be real. Situations with complex <fi can appear as well; then 
the action can contain two non-linear terms, (|0 a | 2 ) 2 and |0 2 | 2 . We will see in Sec. Iri.2.11 
that the phase of a complex order parameter field usually corresponds to a sliding degree 
of freedom of a density wave. 

At positive values of r, the model (fT7|) is in a disordered phase, characterized by 
(0) = 0; at T = the propagator of the fluctuations has a gap of size A(T = 0) = r 
(to zeroth order in go). With decreasing r, this gap decreases, vanishing at the quantum 
phase transition to the ordered phase. In the ordered phase, the magnitude of the order 
parameter is determined by go, at T = at the mean- field level it is (0) = — r / (2g ). For 
N >2 the ordered phase breaks a continuous symmetry and supports gapless Goldstone 
modes with linear dispersion. 

Near the critical point, for N > 2 both amplitude and phase fluctuations can be 
important. However, as will be discussed in Sees. I2.6l and im below the so-called upper- 
critical dimension the critical behavior is determined by phase fluctuations, and in this 
case a theory which neglects amplitude fluctuation - a "hard-spin" theory in contrast 
to the "soft-spin" formulation of the 4 model - is adequate. The hard-spin theory is 
formulated in terms of a field n(x, r) with a unit-length constraint and is known as the 
non-linear sigma model: 



For dimensions 1 < d < 3, equivalently 2 < D < 4 with D — d + z, the 4 model 
()17|) and the non-linear sigma model ()18|) represent different approaches to describe the 
same quantum critical point, and their critical properties are identical [20] ■ Technically, 
the 4 model can be analyzed in an expansion about the upper-critical dimension, i.e., 
for D = 4 — e, whereas the non-linear sigma model allows for an expansion about the 
lower-critical dimension D = 2 + e - those will be briefly discussed in Sec. 12.71 
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Both models (J 17)) and (|18|) possess propagating modes which are undamped at low 
energies, and the quantum transition is characterized by z — 1. In systems with gapless 
fermionic excitations, order parameter fluctuations can, however, decay into particle- 
hole pairs, and this damping requires modifications of the described theories (or even 
different approaches), see Sec. Kill 

As usual, theories of the type (JT7J), supplemented by higher-order terms and/or 
suitable coupling to external fields, can also describe first-order transitions. There, 
the transition point is characterized by the existence of two distinct degenerate global 
minima in the free-energy landscape, which usually leads to a jump in the order 
parameter and other observables when tuning the system through the transition. 
However, no universality is present at first-order transitions, and a continuum 
description is only appropriate for transitions being weakly first order. 

2.5.2. Heisenberg spins. Deriving a continuum model for a system of localized 
quantum-mechanical spins can be conveniently done using spin-coherent state path 
integrals, and the reader is referred to Ref. jH] for an introduction. The orientation 
of each spin is described by a three-component unit vector N, and coherent states |N) 
obey (N|5|N) = 5N. The dynamics of a single spin takes the form of a Berry phase 
term, Sb = J <ir(N|<9 r |N), which is purely imaginary. This term can be shown to be 
equal to iS times the solid angle subtended by the spin vector on the unit sphere (note 
N(/3) = N(0)) - importantly this Berry phase term contains a first-order time derivative. 

The order parameter of interest is again a magnetization density, e.g., for 
antiferromagnets a staggered magnetization. Starting from a lattice spin model, a 
continuum action can be derived which contains gradients of the order parameter field, 
e.g., of the form ([TTj) or (|TSj) . plus a term arising from the spin Berry phases j^j. The 
behavior of resulting low-energy theory crucially depends on those Berry phases, and 
is completely different for ferromagnets and antiferromagnets: In ferromagnetic spin 
systems the Berry phase terms add up and determine the dynamics of the system, as they 
dominate over the second-order time derivative in ()17)) or ()18j) - this is related to the fact 
that the order parameter is a conserved density here. In contrast, in antiferromagnets 
the Berry phase contributions oscillate in sign from site to site, and on average simply 
cancel. Therefore, in many cases Berry phases can essentially be neglected, e.g., in clean 
antiferromagnets on high-dimensional regular lattices, and the resulting model is a 4 or 
non-linear sigma model. However, singular Berry phase configurations can play a role. 
This is crucial in one space dimension, where a corresponding topological term appears 
in the action for half-integer spins, which leads to the well-known gapless "critical" phase 
with power law decay of spin correlations for the d = l Heisenberg model. In contrast, 
for integer spins the topological term is absent, and the system is characterized by an 
excitation gap and conventional short-range spin correlations. 

In a number of cases the analysis of the spin model can be carried out in terms of 
auxiliary particles. One example is the Dyson-Maleev representation of spin operators 
in terms of bosons; another one is the bond-boson representation for spin pairs (dimers) 
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in quantum paramagnets [3]. In both situations one arrives at theories similar to Bose 
gases which will be discussed below. 

2.5.3. Boson and fermion models. Microscopic models of bosons and fermions can have 
a plethora of ordered phases. In many cases, the order parameter can be represented as 
a boson or fermion bilinear, like for spin or charge density waves, and the critical theory 
is of 4 -type with proper modifications (Sec. I3.2.3|) . 

In addition, models of bosons or fermions with a conserved charge Q can have a 
transition between a phase where (Q) is pinned to a quantized value, and a phase where 
Q varies smoothly as function of external parameters. In those cases, the density is 
an appropriate order parameter; however, the theory has to be formulated in terms 
of Bose or Fermi fields and will contain imaginary Berry phase terms for the particle 
dynamics. For instance, varying the chemical potential through zero drives a quantum 
phase transition where the zero-temperature particle density is zero in one phase and 
finite in the other. For bosons in d > 2 this dilute gas quantum critical point is the 
endpoint of a line of finite-T transitions, associated with Bose-Einstein condensation [3J. 

Interesting physics obtains in the presence of strong local interactions, the Hubbard 
model being a paradigmatic example. The physics of the fermionic Hubbard model is 
a subject of vast active research, and cannot be reviewed here (see Sec. 13.41 for some 
remarks on the Mott-Hubbard metal-insulator transition). The bosonic Hubbard model 
is somewhat simpler; it was originally introduced to describe spinless bosons on a lattice, 
representing Cooper pairs undergoing Josephson tunneling between superconducting 
islands, with a strong on-site Coulomb repulsion. This model shows superfluid and 
insulating phases, and an appropriate order parameter is the expectation value of a 
bosonic field, representing the superfluid density (see Sec. 13. 5|) . 

Microscopic realizations of hardcore boson models also arise from models of 
quantum spins, where the degrees of freedom of spins or spin clusters can be re-written 
in terms of bosons with infinite on-site repulsion. An nice example of a dilute Bose 
gas critical point described above is found in quantum antiferromagnets in an external 
field, where the zero-field ground state is a gapped singlet state: Increasing field splits 
the triplet excitation energies, and at a critical field H c the lowest triplet condenses, 
leading to a state with transverse long-range order. Both the T = transition and the 
finite-T transition (present in d > 2), where the transverse order is established, can 
be understood as "Bose-Einstein condensation of magnons". Interestingly, amplitude 
and phase of the condensate wavefunction correspond to magnitude and direction of the 
transverse magnetization, and can thus be measured directly. 

2.6. Fluctuations and critical dimensions 

The critical behavior at a particular transition is crucially determined by the role played 
by fluctuations around the mean- field solution, and their mutual interactions, e.g., the 
go0 4 term in the 4 theory. It turns out that fluctuations become increasingly important 
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if the space dimensionality of the system is reduced. Above a certain dimension, called 
the upper-critical dimension d£, fluctuations are irrelevant, and the critical behavior 
is identical to that predicted by mean- field theory (Ginzburg criterion). For classical 
systems with short-range interactions and a scalar or vector order parameter one finds 
d£ = 4; the mean-field exponents of the corresponding (j) A theory are: 

a = , = 1/2 , 7 = 1, 5 = 3, u = l/2, 77 = . (19) 

Between d~£ and a second, smaller special dimension, called the lower-critical dimension 
d~ , a phase transition exists, but the critical behavior is different from mean- field theory. 
Below d~ fluctuations become so strong that they completely suppress the ordered phase; 
we have d~ = 2 for the classical 4 model with N > 1 components. 

The quantum-classical mapping discussed above implies that for a quantum phase 
transition the critical dimensions are reduced by z compared to the corresponding 
classical transition. Thus, in d = 3 many quantum transitions show mean- field behavior 
(supplemented by logarithmic corrections if z = 1). 

Technically, the existence of the upper-critical dimension is connected with the 
role played by the ultraviolet cutoff of the quantum field theory. For d < d£ the 
integrals arising in perturbation theory are ultraviolet convergent, implying that the 
cutoff can be sent to infinity. Thus, the critical behavior becomes cutoff-independent, 
i.e., truly universal, and observables will follow "naive" scaling, including "tu/T" scaling 
in dynamical quantities, as exploited in Sec. 12 .3\ and exponents will fulfill hyperscaling 
relations (0). Physically, this universality arises from the fact that the non-linear 
interaction terms, e.g., g in the 4 theory, flow towards universal values in the low- 
energy limit. In contrast, for d > d£, quantities will depend on the ultraviolet cutoff, 
and non-universal corrections to observables arise from the non-linearities of the original 
theory - those are termed dangerously irrelevant. Critical exponents will be independent 
of d and thus hyperscaling is violated. 

Importantly, at the critical point in d < d^ the integrals arising in perturbation 
theory are infrared divergent. In other words, (bare) perturbation theory itself is 
divergent below d£ at T = 0. The trick to overcome this problem is the e expansion 
described in the next section. 

For transitions which cannot be easily described by a local order parameter field 
theory, the situation is less clear. For instance, for the disorder- driven metal-insulator 
transition (Sec. 13 .4|) the lower-critical dimension is d~ = 2; it is likely that an upper- 
critical dimension does not exist |2*T] . 

2. 7. Renormalization group approach and calculation of observables 

A powerful tool for the analysis of order-parameter field theories is the renormalization 
group (RG) approach 15 , which is designed to determine the asymptotic low-energy, 
long-wavelength behavior of a system. The momentum-shell formulation of the RG 
proceeds by successively eliminating high-energy degrees of freedom, and at the same 
time renormalizing the couplings of the theory to keep physical properties invariant. 
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Doing so, the coupling constants flow as a function of the cutoff energy, which is 
described by differential RG equations. In the low-energy limit, the RG flow reaches 
a fixed point, and the corresponding fixed point values of the couplings can be trivial 
(zero or infinity) or non-trivial. The trivial values describe stable phases of the system, 
whereas non-trivial finite values correspond to so-called intermediate-coupling fixed 
points, which usually are critical fixed points describing continuous phase transitions. 

Analytical RG computations require a perturbative treatment of the interaction 
terms of the field theory. For d < d£ these grow under RG transformations, i.e., 
interactions are relevant in the RG sense below cf+. For the 4 theory ()17j) the critical 
fixed point with non-linearity g = (the so-called Gaussian fixed point) is unstable w.r.t. 
finite g, i.e., the critical fixed point for d < d£ will be characterized by a finite value 
of g (Wilson-Fisher fixed point). The relevance of the interaction implies that bare 
perturbation theory is divergent at criticality. This problem can be overcome using 
the so-called e expansion: A theory of the (/> 4 -type can be analyzed in an expansion 
around the upper-critical dimension, i.e., in e = c?+ — d. This approach is based on the 
observation that the fixed point value of renormalized coupling g is small near d = d+, 
and a double expansion in g and e allows for controlled calculations. The RG differential 
equations are used to determine RG fixed points and the corresponding values of the 
couplings. For example, in the 4 theory at criticality, the RG equation, describing the 
flow of the dimensionless non-linear coupling g upon reducing the cutoff A, reads 



with a fixed point at g* = (6e)/(N + 8), where e = d+ — d. Physicswise, the finite 
value of g at criticality implies strong self-interactions of the order parameter bosons for 
d < d+; furthermore amplitude fluctuations of the order parameter in a N > 2 4 theory 
are frozen out, and the critical dynamics is carried by phase fluctuations. A different 
expansion can be applied to the non-linear sigma model (j!8j) . namely an expansion 
around the lower-critical dimension d~, i.e., the expansion parameter is e = d — d~ . 
Here, the underlying idea is that the transition temperature t is small near d = d~ , 
and the double expansion is done in t and e around the ordered state. We note that e 
expansions are strictly asymptotic expansions, i.e., the limit e — > has to be taken first, 
and convergence is not guaranteed for any finite e. 

Observables at and near criticality require different treatments depending on 
whether d > d+ or d < d+. Above the upper-critical dimension, bare perturbation 
is usually sufficient. Below d+, one can employ a renormalized perturbation expansion: 
perturbation theory is formulated in terms of renormalized quantities, in the final 
expressions the couplings are replaced by their fixed point values, and the results 
are interpreted as arising from a e expansion of the expected power- law behavior, i.e., 
power laws are obtained by re-exponentiating the perturbation series. The renormalized 
perturbation expansion can thus be understood as a certain resummation technique of 
bare perturbation theory. In some cases, a 1/N expansion, where N is the number 
of order parameter components, can also be used to calculate critical properties for 




(20) 
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arbitrary d. 

The calculation of dynamic and transport properties in the quantum critical regime 
at finite temperatures turns out to be significantly more complicated. Below cf+, 
dynamic properties are dominated by relaxation processes at frequencies u <C T, and 
in this regime the e expansion fails. However, considerable progress can be made 
by applying quasiclassical techniques, where the excitations above the ground state 
can be described by quasiclassical waves or particles [3]. For transport calculations, 
quasiclassical methods are not sufficient, and a quantum-mechanical treatment using, 
e.g., a full quantum Boltzmann equation is necessary. These techniques have been 
applied to a number of models [3], but many issues including the influence of disorder 
remain subject for future work. In contrast, above rf+ perturbation theory is useful 
in many cases, and the presence of quasiparticles allows to analyze transport using a 
quasiclassical Boltzmann equation. 

Beyond perturbative approaches, various numerical techniques have been employed 
to investigate quantum phase transitions, some of which closely follow the general idea 
of renormalization. One example is the numerical renormalization group method (NRG) 
for quantum impurity problems, proposed by Wilson [22]- 

2.8. Example: The transverse-field Ising model 

In this section we want to illustrate the general ideas presented above by discussing the 
Ising model in a transverse field, which displays a paradigmatic example of a quantum 
phase transition. An experimental realization of this model can be found in the low- 
temperature magnetic properties of LiHoF4. This material is an ionic crystal, and at 
sufficiently low temperatures the only magnetic degrees of freedom are the spins of 
the Holmium atoms. They have an easy axis, i.e., they prefer to point up or down 
with respect to a certain crystal axis. Therefore they can be represented by Ising spin 
variables. Spins at different Holmium atoms interact via a dipolar interaction. Without 
external magnetic field the ground state is a fully polarized ferromagnet. 

In 1996 Bitko, Rosenbaum and Aeppli [5j measured the magnetic properties 
of LiHoF4 as a function of temperature and a magnetic field which was applied 
perpendicular to the preferred spin orientation. The resulting phase diagram is shown 
in Fig. |21 A minimal microscopic model for the relevant magnetic degrees of freedom is 
the Ising model in a transverse field. Choosing the z-axis to be the spin easy axis, the 
Hamiltonian is given by 

H = -j'£s;s;-h y E l s? . (21) 

(v) * 

Here Sf and Sf are the z and x components of the Holmium spin at lattice site i, 
respectively. The first term describes the ferromagnetic interaction between the spins, 
while the second term is the transverse magnetic field. For zero field h = 0, the model 
reduces to the well-known classical Ising model, with a ferromagnetic ground state. 
With increasing temperature thermal fluctuations of the spins will reduce the total 
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Figure 2. Magnetic phase diagram of L1H0F4 (after 5 ). The ordered phase breaks 
a Zi symmetry of the Hamiltonian and is bounded by a line of finite-temperature 
phase transitions. This line terminates in the quantum critical point (QCP), where 
ferromagnetic order is destroyed solely by quantum fluctuations. 



magnetization, which will eventually vanish at the critical temperature (about 1.5 K for 
LiHoF 4 ) above which the system becomes paramagnetic. The resulting transition is a 
continuous classical phase transition caused by thermal fluctuations. 

Let us now consider the influence of the transverse magnetic field. Re-writing 
hSf = h/2(Sf + S^), with and being spin flip operators at site i, it is clear that 
the transverse field will cause spin flips which are precisely the quantum fluctuations 
discussed above. If the transverse field becomes larger than some critical field h c (about 
50 kOe in L1H0F4) they will destroy the ferromagnetic long-range order in the system 
even at zero temperature. This transition is a continuous quantum phase transition 
driven exclusively by quantum fluctuations, the phase diagram is of the class shown in 
Fig. Eh. 

For the transverse-field Ising model the quantum-classical mapping discussed in 
Sec. 12.31 can be easily demonstrated at a microscopic level. Consider a one-dimensional 
classical Ising chain with the Hamiltonian 

N 

H d = -Jj2S*S? +1 . (22) 

i=i 

Its partition function is given by 

N N 

Z = Tr e~ H ^ T = Tr exp(J £ S*S* +1 /T) = Tr ]J M i4+1 , (23) 

i=i i=i 

where M is the so called transfer matrix. It can be represented as 

/ JIT — JIT \ 

M = 6 _ J/T 6 J/T = e J ' T {\ + e~ 2J l T S x ) « exp(J/T + e~ 2J l T S x ) . 
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Except for a multiplicative constant the partition function of the classical Ising chain 
has the same form as that of a single quantum spin in a transverse field, Hq = —h S x , 
which can be written as 

N 

Z = Tt e- H ^ = Tr e hSX ^ =Trf[ e hSX '^ . (24) 

i=l 

Thus, a single quantum spin can be mapped onto a classical Ising chain. Under this 
mapping, the classical temperature is related to a coupling constant in the quantum 
model, h = e~ 2J ^ T (being the inverse correlation length of the classical problem) , and the 
classical system size, N, maps onto the inverse of the quantum temperature, Tq. These 
considerations can easily be generalized to a <i-dimensional transverse-field (quantum) 
Ising model which can be mapped onto a (d+ l)-dimensional classical Ising model. 
Consequently, the dynamic exponent is z = 1 here. 

3. Bulk quantum phase transitions 

In this section we will review a number of classes of quantum phase transitions occuring 
in condensed matter systems. We will focus our attention on continuous quantum phase 
transitions, which are characterized by diverging correlation lengths and times. 

The continuous transitions discussed here are so-called bulk transitions, meaning 
that degrees of freedom of the whole sample become critical at the transition point, as is 
the case with most familiar phase transitions. This implies that the singular part of free 
energy of the system scales as L d where L is the linear spatial extension of the system 
and d the space dimensionality. In Sec. |Uwe will discuss so-called boundary transitions 
where only a subset of the available degrees of freedom becomes critical. 

Most of the following considerations will deal with "clean" systems, where 
translational invariance is unbroken and quenched disorder is absent. The interplay 
between disorder and criticality is a notoriously difficult issue - even in classical systems 
- and only understood in some specific models to date; we will touch upon this issue in 
Sec. EH 

3.1. Quantum phase transitions and fermions 

Quantum critical behavior depends crucially on whether order parameter fluctuations 
can couple to low-energy fermionic excitations (or, more generally, to other non-critical 
soft modes of the system). In the absence of gapless fermions the order parameter 
fluctuations are the only low-energy excitations in the vicinity of the critical point, and 
the critical theory can be formulated in terms of the order parameter alone. Technically, 
the step of integrating out the fermions from the action describing the interacting 
electron system does not lead to any divergencies. 

Critical dynamics can be fundamentally changed if order parameter fluctuations 
couple to low-energy fermions, as are present in metals. Here, integrating out fermions 
can lead to divergencies. In some cases, it is believed that a formulation of a critical 
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theory for the order parameter alone is still possible, however, the presence of low- 
energy fermions leads to non-analytic terms in the effective order parameter theory - 
an example is the Landau damping of critical magnetic fluctuations in metals. 

In the following, we specialize the discussion to conventional order parameters which 
can be written in terms of local fermion operators and carry a certain momentum 
Q. In a translational invariant system, momentum conservation dictates that an 
effective coupling between order parameter and fermionic excitations is only possible 
if the characteristic wavevector Q of the order parameter connects Fermi surface points, 
as fermions can be scattered off order parameter fluctuations acquiring a momentum 
change of Q. Four cases have to be distinguished, characterizing the order parameter 
momentum Q: (i) If the order parameter has zero total momentum, as for a ferromagnet, 
then it couples to particle-hole pairs at the whole Fermi surface. We will further discuss 
this case in Sec. 13.2.31 For an order parameter with finite momentum Q, the remaining 
cases are: (ii) Q cannot connect points on the Fermi surface. Then low-energy fermions 
are scattered into states with higher energy, perturbation theory in the corresponding 
coupling is convergent and only gives a renormalization of parameters. Thus, the critical 
behavior is the same as in insulators, (iii) Q connects a (d— 2)-dimensional set of points, 
i.e., lines in d— 3 or discrete points in d— 2. This leads to damping of the order parameter 
modes due to the decay into particle- hole pairs. At least in d = 3 it is believed that this 
effect can be fully captured by a Landau damping term in the effective action, which 
increases the dynamic critical exponent; explicit examples are in Sec. 13.2.31 (iv) For 
so-called perfect nesting, Q can connect whole parts of the Fermi surface, and then a 
transition generically occurs at infinitely small coupling. 

An effective coupling between order parameter fluctuations and low-energy fermions 
can also be present in unconventional superconductors which have point or line nodes 
in their gap function, because then a Fermi surface still exists, although it is no longer 
{d—1) -dimensional. A particular example is a two-dimensional <i-wave superconductor, 
see Sec. 13.2.41 Similar considerations apply to semimetals. 

For unconventional order parameters, e.g., for glassy order or topological order, the 
situation is less clear. It has been proposed that fluctuations of those order parameters 
can couple to fermions on the entire Brillouin zone, leading to various scenarios of "local 
criticality" . 

In general, a critical point involving both low-energy order parameter and fermionic 
modes can be expected to display two time scales, because order parameter fluctuations 
can relax both via (damped) order parameter modes and via fermions. Consequently, 
the critical theory should be characterized by two dynamic critical exponents |23j . A 
proper RG analysis of such a theory, keeping both low-energy degrees of freedom, has 
so far only been carried out for the ferromagnetic transition |24j . 
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3.2. Phase transitions with conventional order parameters 

This section describes the physics of continuous quantum phase transitions with local 
order parameters, involving, e.g., spin, charge, Peierls, or phase order. In space 
dimensions larger than the classical lower-critical dimension (see Sec. 12. 6|) . those 
transitions are zero-temperature endpoints of a line of classical phase transitions, Fig.^D. 

3.2.1. Order parameters. Many order parameters can be expressed as expectation value 
of combinations of particle creation and annihilation operators which are local in space 
and time. The correlation functions of the so-defined local operator 0(Rj,£) do not 
decay to zero in the long-distance limit in the ordered phase (long-range order). As 
the local quantity can oscillate in space, like the magnetization in an antiferromagnet, 
the order parameter 0(R») is usually defined as the expectation value of O with the 
microscopic oscillations removed, 



and Q is the ordering wavevector or characteristic momentum. Sufficiently close to the 
ordering transition, 0(Rj) is now slowly varying in space, and can be replaced by a 
continuum field 0(r). In a path integral description, we then arrive at a bosonic field (f) 
carrying momentum Q. 

The order parameter can have multiple components, i.e., be a spinor, vector, or 
tensor, and the underlying microscopies then defines the symmetry properties. The 
symmetry in turn determines the allowed terms in the low-energy action and thus the 
universality class. It also decides about the existence of low-energy modes: an ordered 
phase which breaks a continuous symmetry supports Goldstone modes. In addition, 
the order parameter can carry charge, e.g., if it describes particle-particle rather than 
particle-hole pairing; the presence of long-range Coulomb interactions can then change 
the nature of both the collective modes and the transition. 

We list a few common examples for order parameters below. For density waves, the 
number of components is N = 1, and the density oscillates around its average value as 



Similar scalar or Ising order parameters also occur, e.g., for alloy ordering and for 
Peierls transitions - in the latter the system undergoes a spontaneous doubling of the 
unit cell, induced for instance by strong electron-phonon or spin-spin interactions. For 
spin density waves, we have a vector order parameter, N = 3 and a — x,y, z, and the 
spin density modulation is given by 



Here Q s = describes a ferromagnet. In the presence of strong spin anisotropies, 
fluctuations in some directions can be frozen out at low energies, then leading to 
magnetic theories with N = 1 or iV = 2. XY rotor models (N = 2) can also represent 
phase variables, appropriate, e.g., for phase fluctuations in superconductors. 




(25) 




(26) 




(27) 
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Interestingly, in Eqs. (12612 7|) is real only for the cases Q = and Q = (tt, 7T, . . .) 
(assuming a hypercubic lattice with unit spacing), otherwise <fi can be complex with 
the phase describing a sliding degree of freedom of the density wave; this phase takes 
discrete values for commensurate wavevectors. In the magnetic case, both collinear and 
non-collinear spin density waves can be described by Eq. (|2*Tj) [25 : 

collinear: <p sa = e t0 n a with n a real , 

spiral: <p sa = rii a + m 2Q with ri\^a real and ni a n 2a = . (28) 

We remind the reader that in spin systems Berry phases are important, which can render 
the analysis in terms of a slowly varying order parameter alone invalid. 

Further interesting order parameters are those with an unconventional symmetry; 
for fermions bilinears this implies that different parts of the Fermi surface contribute 
with different weights or signs. Examples are superconductors with p, d, or higher 
angular momentum pairing symmetries. Also particle-hole pairing of unconventional 
symmetry is possible; in the charge channel this can be understood as deformation of 
the Fermi surface leading to "nematic" states - the transition of a Fermi liquid into such 
a state is also known as Pomeranchuk instability. While most of the listed cases have 
wavevector Q = 0, finite Q cases are possible as well, an example being the so-called 
d-density wave or staggered flux order |2S1, with Q = (n, tc). 

3.2.2. Insulators. The simplest quantum-phase transitions are those where the order 
parameter fluctuations do not couple to low-energy fermionic excitations - this is 
generically the case in insulators where the gap to charge excitations is large. 

Then, it is sufficient to consider a theory of the order parameter fluctuations alone. 
In most cases, the appropriate continuum theory takes the form of a 4 theory (|17|) 
described in Sec. 12 ,5[ with a dynamic critical exponent z — 1. For clean systems, results 
for static quantities at T = can be obtained from the corresponding classical transition 
via the quantum-classical mapping. Most critical properties of such transitions are well 
understood and can be found in textbooks jTTJ [T2] ; we will not describe them further 
here. In contrast, real-time dynamics and transport, in particular at finite temperatures, 
cannot be easily extracted from the classical results, and are demanding subjects of 
current research [3 . For disordered systems a number of complications arise, described 
in Sec. 13.31 and only relatively few results are available. 

3.2.3. Metals: Non-Fermi liquid behavior. Fermi liquids in space dimensions d > 2 
can display ground-state instabilities towards states with order parameters described in 
Sec. I3.2.1[ which are of considerable practical interest. In most cases, fermionic low- 
energy modes couple to the order parameter fluctuations; many aspects of the resulting 
theories are not fully understood to date. In d > 2 the quantum critical theories are 
then above their upper-critical dimension, therefore non-universal features abound, and 
microscopic details like the Fermi surface topology become important. A unified picture 
of these transitions is not yet available. 
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The investigation of quantum phase transitions in Fermi liquids was pioneered by 
Hertz [27] who studied the ferromagnetic transition by a renormalization group method. 
Millis j2HJ considerably extended this work and computed finite-temperature crossover 
functions for several magnetic transitions. The theoretical approach (commonly called 
Hertz-Millis theory) proceeds by integrating out the fermions from the full interacting 
theory, keeping only the low-energy order parameter fluctuations in the theory. Although 
this is a formally exact step, the coefficients of the resulting action will in general be 
non-local in space and time because the fermionic excitations are gapless. Assuming 
that the action can be expanded in powers and gradients of the order parameter field 
and assuming that it is permissible to truncate such an expansion, one obtains 



for the case that the ordering wavevector Q connects a (d— 2)-dimensional set of points 
of the Fermi surface. Compared to the usual 4 theory, the dynamical term is now 
replaced by \u n \ describing the Landau damping of the order parameter fluctuations, 
i.e., the possibility to decay into particle-hole pairs. The theory ()29|) has dynamic 
exponent z = 2; for the case of a ferromagnet (Q = 0) the damping term gets replaced 
by \ui n /k\, hence z = 3. Note that fermionic gap, which occurs in the ordered phase 
near the momenta connected by Q, is not captured by (|29)L therefore this action is not 
to be taken seriously at r < 0, T = [Sj. 

The type of quantum critical behavior within Hertz-Millis theory [2H] is dictated 
by the fact that the interaction term is formally irrelevant above the upper-critical 
dimension, therefore only singular corrections to a Gaussian theory emerge. For 
example, the specific heat of the three-dimensional antiferromagnet behaves as C /T = 
7 + A\fT in the quantum critical region, in other words, the specific heat coefficient 7 
of the Fermi liquid does not diverge upon approaching the quantum phase transition. 

Subsequent theoretical work has indicated that the above assumption of a regular 
gradient expansion of the action breaks down in a number of important cases. In 
the antiferromagnet case, an analysis of the spin-fermion model has indicated that the 
Hertz-Millis theory is invalid in d—2 dimensions [29 . For the ferromagnetic transition, 
it has been shown that non-analytic terms in the spin susceptibility of a usual Fermi 
liquid lead to corresponding corrections in the effective order parameter theory [19J, 
which drastically modify the critical behavior. Notably, the resulting interaction terms 
are non-local, and the explicit form of the four-point vertex is not known, rendering 
explicit calculations very difficult in this approach. Clearly, an improved theory should 
be local and should explicitely keep both the order parameter fluctuations and the low- 
energy fermionic modes, making manifest the presence of two time scales in the effective 
theory. Such a calculation has been recently performed for the ferromagnetic transition 
in metals. (In an unconventional superconductor, the situation is somewhat simpler 
due to the absence of a full Fermi surface, and both order parameter and fermionic 




(29) 
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modes can be kept, see Sec. 13.2.41 ) We briefly describe the results for the ferromagnetic 
transition [21]: In a clean metal with 1 < d < 3 the singular behavior of the Fermi 
liquid destabilizes any continuous transition with z = 3, therefore only a first-order 
transition to a ferromagnetic state or a z = 2 transition to a state with spiral order 
appear possible. In d = 3 the conclusions are similar, however, the corrections to the 
z = 3 critical point are only logarithmic, and in certain parameter regimes a continuous 
transition to a ferromagnetic state may be restored. In the disordered case, i.e., for 
diffusive spin dynamics, it was found that the z = 4 fixed point of Hertz is unstable, and 
instead a new transition with logarithmic corrections to a Gaussian fixed point with 
z — d occurs. 

On the application side, itinerant ferromagnetic transitions have been found in 
MnSi and other intermetallic compounds; in many cases those transitions become 
first order at low temperatures, which appears consistent with the theoretical 
findings. Antiferromagnetic transitions appear in both transition-metal and rare-earth 
compounds. Great interest has been focussed on the normal-state properties of the 
high-temperature superconductors, where simple square-lattice antiferromagnetism as 
well as collinear spin and associated charge density waves ( "stripes" ) have been detected 
[001 ED E21 EH] • Those order parameters apparently coexist and/or compete with <i-wave 
superconductivity at low temperatures. Also, phase transitions involving circulating 
currents have been proposed to explain the unusual normal-state properties of high-T c 
cuprates [201 El] • 

Much experimental and theoretical work has been performed on heavy-fermion 
compounds. A well-studied compound is CeCug.^Au^, which displays a transition 
between a heavy Fermi liquid phase and an antiferromagnetic metallic phase at x — 0.1. 
Interestingly, at the transition the specific heat coefficient diverges as C/T oc In T, and 
neutron scattering indicates uj/T scaling in the response functions [00]. However, such 
scaling is only expected below the upper-critical dimension, see Sec. I2.6[ and is clearly 
at odds with the Hertz-Millis approach. Some proposals have been made to resolve 
this inconsistency, relating the observed quantum critical behavior to the breakdown 
of Kondo screening [001 EH EH] • In Ref. |0Z], a theory of "local" critical behavior has 
been proposed in the spirit of the dynamical mean- field approach |39j : however, non- 
trivial behavior emerges here only for two-dimensional systems. A different approach 
[OH] has associated the suppression of Kondo screening with the emergence of a non- 
magnetic "fractionalized Fermi liquid" phase, obtained by suppressing magnetic order 
in the regime of weak Kondo screening; see also Sec. 13.6.21 

3.2.4- Superconductors. It is conceivable that a quantum phase transition occurs with 
superconductivity being present on both sides including the transition point. The main 
difference to the metallic case discussed above is that low-energy fermions are gapped 
due to the presence of the superconducting gap. Then, in general the order parameters 
fluctuations are undamped, and the phase transitions are in the same universality class 
as the ones in insulators. 
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However, in unconventional superconductors showing gapless points or lines in 
momentum space (called nodal points in the following), damping can still occur. An 
effective coupling between order parameter fluctuation and fermions is only present if 
the wavevector Q of the order parameter fluctuations does connect two nodal points, 
see Sec. 13.11 Those special cases are ideally treated in a theory which keeps both 
order parameter and fermionic excitations, and examples have been considered in Refs. 
[4*0*1 "4*T* . Candidate transitions are those with a particular (fine-tuned) finite wavevector 
Q, like spin or charge density waves, and transitions with Q = which do not require any 
fine tuning. Examples in the latter class are particle-hole pairing transitions which lead 
to "nematic" states, and transitions involving the onset of secondary superconducting 
pairing |41j . Note that for those transitions the onset of the secondary pairing occurs 
at a finite value of the corresponding coupling because the density of states in the 
background superconducting state vanishes at the Fermi level. 

We briefly sketch the theory of the transition between two superconducting states 
with d x 2_ y 2 and d x 2_ y 2 + id xy pairing symmetry [IT] - this is one of the simplest models 
where both order parameter fluctuations and low-energy fermionic excitations can be 
treated on equal footing. Importantly, the fluctuations of the secondary d xy order 
parameter are represented by an Ising field (0), because the fluctuations of its phase 
relative to the background <i-wave pairing are massive. Thus we simply have 



Sa, = Id xdr 



kd T <pT + c ^(v x <p) 2 + r -<p 2 + 9 -Z~ 



(30) 



2 V r/ 2 v 2 r 4! 

as in Eq. ([17)1 . We explicitely include the low-energy fermions which have a linear 
spectrum; for one pair of nodes their action can be written as 

r d 2 k 

Sm = Jj^TJ2 *Ia (~«"n + v F k x r z + v A k y r x ) * la . (31) 

Here ^>i a = (fi a , Sabflb) is a Nambu spinor of fermions from opposite nodes (1 and 3) of 
the ci-wave superconductor, r a are Pauli matrices which act in the fermionic particle- 
hole space, k XtV measure the wavevector from the nodal points and have been rotated 
by 45 degrees from q X)V co-ordinates, and vp, va are velocities. (An analogous term 
describes the nodes 2 and 4.) The full action is then 5 = 5^ + S^i + S^ 2 + S^^ where 
the final term in the action, S^, couples the bosonic and fermionic degrees of freedom: 

= Jd 2 xdr[X <P (V la Mi^i a + ^ 2a M^ 2a ) ] , (32) 

where M 12 are now specific for the secondary order parameter considered, for d xy pairing 
we have Mi = r y , M 2 = —r y |UJ. This theory can be analyzed by perturbative RG 
together with e expansion in standard fashion. The resulting transition has dynamical 
exponent z — 1, and the ordered phase breaks a Z 2 symmetry. However, the transition 
is not in the Ising universality class due to the coupling to the fermions. An explicit 
calculation of the boson propagator shows that it contains non-analytic terms in both 
k and u which are different from simple Landau damping. The sketched theory has 
been proposed to describe the strong damping of nodal fermions in cuprates below the 
superconducting T c [4"T] . 



QUANTUM PHASE TRANSITIONS 



28 



3. 3. Quantum phase transitions and disorder 

In application to real system the influence of static or quenched disorder on the 
properties of a quantum phase transition is an important aspect. Remarkably, the 
effect of disorder is not completely understood even for classical phase transitions. 

In a theoretical description, quenched disorder can occur in different ways: on a 
microscopic level, e.g., random site energies or bond couplings, or randomly distributed 
scattering centers are possible. In an order parameter field theory, disorder usually 
translates into a random mass term for the order parameter fluctuations. Importantly, 
applying the quantum-classical mapping (Sec. 12.3)) to a quantum problem with quenched 
disorder leads to a (c?+z)-dimensional field theory with strongly anisotropic correlated 
disorder because the disorder is frozen in the time direction. In some cases lattice 
effects not captured by the field theory can be important, this applies, e.g., to all types 
of percolation problems. Moreover, disordering a quantum model can lead to random 
Berry phase terms which have no classical analogue, an example are diluted Heisenberg 
ant if erromagnet s . 

A basis for the discussion of the effect of disorder on continuous transitions is given 
by the Harris criterion ^21 S3] , which states that if the correlation length exponent v of 
a given phase transition obeys the inequality v > 2/d, with d the space dimensionality 
of the system, then the critical behavior is unaffected by quenched disorder. In the 
opposite case, v < 2/d, the disorder modifies the critical behavior, leading either (i) 
to a new critical point that has a correlation length exponent v > 2/d and is thus 
stable, or (ii) to an unconventional critical point where the usual classification in terms 
of power-law critical exponents looses its meaning, or (iii) to the destruction of a sharp 
phase transition. The first possibility is realized in the conventional theory of random- 
ly classical ferromagnets [44], and the second one is probably realized in classical 
ferromagnets in a random field jIH]- The third one has occasionally been attributed 
to the exactly solved McCoy- Wu model |%B] . This is misleading, however, as has been 
emphasized in Ref. jUJj; there is a sharp, albeit unorthodox, transition in that model, 
and it thus belongs to category (ii). 

For first-order transitions, Imry and Wortis [47] have proposed a criterion stating 
that below a certain critical dimension (which is usually 2) the transition will be 
significantly rounded due to domain formation. The question of whether the possibly 
resulting continuous transition shows universal scaling is not settled [H]. 

Independent of the question of if and how the critical behavior is affected, disorder 
leads to very interesting phenomena as a phase transition is approached. Disorder in 
classical systems generically decreases the critical temperature T c from its clean value 
T c °. In the temperature region T c < T < T c ° the system does not display global order, but 
in an infinite system one will find arbitrarily large regions that are devoid of impurities, 
and hence show local order, with a small but non-zero probability that usually decreases 
exponentially with the size of the region. These static disorder fluctuations are known 
as "rare regions", and the order parameter fluctuations induced by them as "local 



QUANTUM PHASE TRANSITIONS 



29 



moments" or "instantons" . Since they are weakly coupled, and flipping them requires 
to change the order parameter in a whole region, the local moments have very slow 
dynamics. Griffiths [SD] was the first to show that they lead to a non-analytic free 
energy everywhere in the region T c < T < T c °, which is known as the Griffiths phase, 
or, more appropriately, the Griffiths region. In generic classical systems this is a weak 
effect, since the singularity in the free energy is only an essential one. (An important 
exception is the McCoy- Wu model [35].) 

Turning now to quantum phase transitions, it is expected that disorder has a 
stronger effect compared to the classical case due to the disorder correlations in time 
direction. Consequently, the Harris criterion is still given by v < 2/d (not d + z). 
A prototypical and well-studied model is the quantum Ising chain in a transverse 
random field, investigated by Fisher jlHj and others, where the phase transition is 
controlled by a so-called infinite-randomness fixed point, characterized by a wide 
distribution of couplings and activated rather than power-law critical behavior. The 
Griffiths singularities are also enhanced compared to the classical case: a number of 
observables display power-law singularities with continuously varying exponents over a 
finite parameter region in the disordered phase. Numerical simulations [53 j suggest 
that random-singlet phases and activated criticality may not be restricted to one- 
dimensional systems, raising the possibility that exotic critical behavior dominated 
by rare regions may be common to certain quenched-disorder quantum systems, in 
particular those with Ising symmetry. Recent investigations of two-dimensional diluted 
bilayer antiferromagnets with Heisenberg symmetry indicate conventional critical 
behavior, and in addition an interesting interplay of quantum and geometric criticality 
at the percolation threshold |5T] . 

All the examples mentioned so far display undamped order parameter dynamics 
(z — 1 in the clean system). In contrast, in itinerant electron systems the dynamics 
is overdamped due the coupling to fermions (z > 1 in the clean limit, see Sees. 13.11 
and 13.2.5)1 . which also changes the behavior upon introduction of disorder. It has been 
shown recently that for Ising order parameter symmetry and overdamped dynamics a 
continuous transition is always rounded due to the presence of ordered islands arising 
from a 1/r 2 interaction in the effective classical model [22] • For continuous symmetries 
it is likely that a sharp transition survives. A prototype example is the disordered 
itinerant antiferromagnet, where it has been shown that a finite-disorder fixed point 
is destabilized by the effects of rare regions. The RG indicates run-away flow to large 
disorder, and definite conclusions regarding the transition cannot be drawn 54 . In 
the context of heavy-fermion compounds, it has been suggested that non-Fermi liquid 
behavior can occur in the quantum Griffiths region of a disordered Kondo lattice model 

Systems which are dominated by disorder and frustration can have a variety of new 
stable phases not known from clean systems, e.g., spin or charge glasses and various 
types of infinite-randomness phases. Associated phase transitions then belong to new 
universality classes, an interesting example being the quantum spin glass transition, i.e., 
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the T = transition between a paramagnet and a spin glass Such transitions may be 
relevant for certain strongly disordered heavy-fermion compounds, but also speculations 
in the context of high-T c superconductivity have been put forward (HE] • 

3.4- Metal-insulator transitions 

Metal-insulator transitions are a particularly fascinating and only incompletely 
understood class of quantum phase transitions. Conceptually, one distinguishes between 
transitions in models of noninteracting electrons, arising from solely from lattice effects, 
and transitions of interacting electrons. Examples in the first class are band or Peierls 
metal-insulator transitions, as well as the disorder- driven Anderson transition. The 
prominent example in the second class is the Mott transition of clean interacting 
electrons. At the Anderson transition, the electronic charge diffusivity D is driven to zero 
by quenched disorder, while the thermodynamic properties do not show critical behavior. 
In contrast, at the Mott transition the thermodynamic density susceptibility dn/dfi 
vanishes due to electron-electron interaction effects. In either case, the conductivity 
a = (dn/dfi)D vanishes at the metal-insulator transition. It is worth emphasizing that 
a sharp distinction between metal and insulator is possible only at T = 0; in some cases, a 
first- order finite-temperature transition between a "good" and a "bad" metal can occur. 

Let us briefly discuss the disorder-driven metal-insulator transition. Introducing 
quenched disorder into a metallic system, e.g., by adding impurity atoms, can change 
the nature of the electronic states from spatially extended to localized jHZj. This 
localization transition of disordered non-interacting electrons, the Anderson transition, 
is comparatively well understood, see Ref. [HE!- The scaling theory of localization [HH] is 
based on the assumption of hyperscaling and predicts that in the absence of spin-orbit 
coupling or magnetic fields all states are localized in one and two space dimensions 
for arbitrarily weak disorder. Thus, no true metallic phase exists for al = 1,2. In 
contrast, in three dimensions there is a phase transition from extended states for weak 
disorder to localized states for strong disorder. These results of the scaling theory are in 
agreement with large-scale computer simulations of non-interacting disordered electrons. 
The effect of weak disorder can be captured by perturbation theory, leading to the well- 
known "weak-localization" quantum corrections to the conductivity, which are divergent 
in d — 1,2. A field-theoretic description of the Anderson transition has been pioneered 
by Wegner 60J. Starting from electrons in a random potential, the disorder average is 
performed using the Replica trick, and the arising interaction terms can be decoupled 
using complex matrix fields. The theory then takes the form of a non-linear sigma model 
in terms of these matrix fields; it has been subsequently analyzed at the saddle-point 
level and by RG methods, with results consistent with the scaling theory [2"T] . 

A different limiting case is the interaction-driven Mott transition of clean electrons 
which can occur at commensurate band fillings. The most studied model showing such 
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a transition is the one-band Hubbard model, with the Hamiltonian 

H = ~ Yl tiMvCja + c )* c io) + U Y1 n ^ n H ' ( 33 ) 
(ij)a i 

defined on a regular lattice, with hopping energies and an on-site Coulomb repulsion 
U. At half filling, the system is metallic at U = 0, but insulating at U —>■ oo because large 
U prevents doubly occupied sites and thus completely suppresses hopping. In general, 
an interaction-driven metal-insulator transition may be accompanied by the occurrence 
of magnetic order (Mott-Heisenberg transition). In contrast, the transition from a 
paramagnetic metal to a paramagnetic insulator is termed Mott-Hubbard transition. 
Non-perturbative techniques are necesssary to analyze the Mott transition because it 
generically occurs at intermediate coupling (except for perfect nesting situations). 

The Hubbard model ()33|). both at and away from half-filling, has been extensively 
studied in the context of high-temperature superconductivity, and a plethora of phases 
have been proposed as function of parameters and doping. However, many properties 
of this apparently simple one-band model in two and three dimensions are still under 
debate, as analytical calculations often have to rely on uncontrolled approximations, 
and numerical methods are limited to very small system sizes. In one dimension, 
special techniques have provided a number of exact answers, see Ref. [HI] . Beyond 
that, considerable progress has been made in the framework of dynamical mean-field 
theory (DMFT) [21], which represents a local approximation to the self-energy and 
becomes exact in the limit of infinite space dimensions. In this approach, the lattice 
model is mapped onto a quantum impurity model supplemented by a self-consistency 
condition. Extensive numerical studies have established that the T = Mott transition 
in the DMFT limit is characterized by the continuous vanishing of the quasiparticle 
weight upon increasing U/t; however, the insulator displays a preformed gap, and no 
dynamic critical behavior is present [23 E2] • 

In reality, both Coulomb interaction and disorder are present, and their interplay is 
poorly understood. The conventional approach to the problem of disordered interacting 
electrons is based on a perturbative treatment of both disorder and interactions [oil] I64j . 
It leads to a scaling theory and a related field-theoretic formulation of the problem [65J, 
which is a generalization of Wegner's theory for non-interacting electrons jHO], and was 
later investigated in great detail within the framework of the renormalization group 
(for a review see Ref. [21]). One of the main results is that in the absence of external 
symmetry breaking (spin-orbit coupling or magnetic impurities, or a magnetic field) a 
phase transition between a normal metal and an insulator only exists in dimensions 
larger than two, as was the case for non-interacting electrons. In two dimensions 
the results of this approach are inconclusive since the renormalization group displays 
runaway flow to zero disorder but infinite interactions. Furthermore, it has not been 
investigated so far, whether effects of rare regions (Sec. 13. 3|) would change the above 
conclusions about the metal-insulator transition. It is known that local moments tend 
to form in the vicinity of the transition, however, their role in the critical properties 
has not been clarified. Conceptually, it is clear that a perturbative approach is not 
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applicable in a regime of strong interactions. 

Experimental work on the disorder-driven metal-insulator transition, mostly on 
doped semiconductors, carried out before 1994 essentially confirmed the existence of a 
transition in three dimensions jHSj while no transition was found in two-dimensional 
systems. Therefore it came as a surprise when experiments on Si-MOSFETs [TU] 
revealed indications of a true metal-insulator transition in two dimensions. These 
experiments have been performed on high-mobility samples with a low carrier density. 
Therefore, electronic interactions are strong: for an electron density of 10 11 cm -2 the 
typical Coulomb energy is about 10 meV while the Fermi energy is only about 0.5 
meV. Interaction effects are a likely reason for this new behavior in two dimensions, 
however, a complete understanding has not yet been obtained. It is not settled whether 
the observations are associated with a true zero-temperature transition, or whether 
they represent a intermediate-temperature crossover phenomenon. Recent experiments 
indicate glassy charge dynamics on both the metallic and the insulating side of the 
apparent transition, further pointing to an interesting interplay between disorder and 
interactions. 

A number of ideas developed in context of the Anderson transition also apply to 
transitions in quantum Hall systems. In each quantum Hall plateau phase, electrons are 
localized due to disorder, and the longitudinal T = conductivity vanishes. Transitions 
between different quantum Hall states involve a diverging correlation length and scaling 
behavior of the conductivity. In addition to scaling theories, approaches based on 
quantum tunneling network models, as first proposed by Chalker and Coddington, have 
been successful in explaining a number of experimental findings. However, at present the 
role of electron interactions at the quantum Hall critical points has not been thoroughly 
clarified. For a full exposure we refer to recent review articles [Tl I67j. 

Summarizing, metal-insulator transitions, in particular those occurring at 
intermediate or strong interactions, are in the focus of current research activity, and 
many aspects, including the interplay with lattice effects and magnetism, remain to be 
understood. 

3.5. Superfluid-insulator transitions 

Transitions between superfluid and insulating states can, like metal-insulator 
transitions, be driven by disorder or by strong Coulomb interaction; in the case of a 
charged superfluid, i.e., a superconductor, applying a magnetic field can also lead to a 
superconductor-insulator transition. Both the disorder and field-driven transitions have 
been extensively studied in the context of thin superconducting films. Interaction-driven 
(Mott) transitions can be realized in Josephson junction arrays where the ratio between 
the inter-island Josephson coupling and the intra-island Coulomb interaction can be 
varied [I]. Very recently, clear signatures of a Mott superfluid-insulator transition have 
been observed in a system of ultracold atoms in an optical lattice (68) . 

Usually, models of bosons are employed to describe the superfluid. For the disorder- 
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free case with interactions one is lead to the boson Hubbard model, introduced in 
Sec. 12.51 which contains a tight-binding hopping energy t and an on-site repulsion 
U. The ground-state phase diagram is well understood: In a plane parametrized by 
the chemical potential fi/t and the ratio t/U the large-t regime is superfluid whereas 
in the small-t regime Mott phases occur where the density is fixed to integer values. 
The transitions are described by an 0(2) rotor model capturing the condensate phase 
dynamics provided the density is fixed - this is the case on the tips of the so-called 
Mott lobes of the phase diagram. Otherwise, the complex condensate density itself is 
the appropriate order parameter, and the transition is in the dilute Bose gas class jSj. 

In the presence of disorder, the situation is somewhat more complicated. In addition 
to the Mott insulator and superfluid phases a Bose glass phase appears, which is believed 
to be insulating. Particular attention has been focussed on the two-dimensional case 
where the scaling dimension of the conductivity turns out to be zero. Thus, scaling 
theories predict a universal metallic resistance exactly at the transition point between 
superconductor and insulator An important feature of the theory is the relevance 

of long-range Coulomb interaction at the transition, which leads to a dynamical critical 
exponent being unity (HS]- A duality description of the two phases, involving bosons and 
vortices, has been put forward: in the superconducting phase the bosons are condensed 
and the vortices localized whereas in the insulating (Bose glass) phase the bosons are 
localized and the vortices condensed - this corresponds to a superfluid of vortices [7H] . 

While these theories were confirmed by a number of experimental studies, they have 
been called into question by recent experiments on superconducting films in a field, which 
suggest the existence of a finite intermediate regime of metallic behavior. A number 
of proposal have been put forward, ranging from a metallic Bose glass, or a true Bose 
metal, to metallicity produced by the influence of dissipation [71]. Other speculations 
concern the possible influence of fermionic excitations, which are usually neglected in 
boson-only models. In this context one has to keep in mind that in two dimensions 
a fermionic metallic phase is believed to be absent due to weak localization; however, 
the experimental findings on two-dimensional electron gases mentioned in Sec. 13.41 have 
called this into question as well. Thus, the fascinating fields of metal-insulator and 
superconductor-insulator transitions are firmly connected, and the interplay between 
disorder and interactions remains an important topic of future research. 

3.6. Phase transitions involving topological order 

Phase transitions can involve the onset of so-called topological order; in this case the 
ordered state is not characterized by a local order parameter; however, the two phases 
can be distinguished from their global properties. In many cases, topological order can 
be rephrased as the suppression of topological defects, either in physical (e.g. spin) 
or emergent (gauge field) degrees of freedom. The transition from a topologically 
ordered to a disordered state can be described as the proliferation of topological defects. 
Topological order occurs in different contexts, e.g., in the quantum Hall effect, in 
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superconductors, and in spin systems; we will describe a few interesting transitions 
below. 

3.6.1. Kosterlitz-Thouless transition. A well-known classical example is the finite- 
temperature transition of the d=2 XY model, i.e., a model where spins are confined 
to the x-y plane which is equivalent to a 0(2) rotor model. At high temperatures 
the spins are disordered with exponentially decaying correlations, whereas at low, non- 
zero temperatures the spin display quasi-long-range order with power-law correlations. 
This quasi-ordered phase has topological vortex excitations which are suppressed at 
low temperatures, i.e., the density of vortex-antivortex pairs is small. The two phases 
are separated by a Kosterlitz-Thouless transition at Tkt where vortex pairs unbind 
upon increasing temperature; the transition is accompanied by a universal jump in 
the superfluid density. The XY model is of relevance for the phase dynamics of 
two-dimensional superconductors and has been discussed extensively in the context of 
underdoped high-temperature superconducting cuprates. 

A Kosterlitz-Thouless quantum phase transition occurs at T = in the (1 + 1)- 
dimensional XY model, or equivalently in the 0(2) rotor model in d— 1 - this follows 
straightforwardly from the quantum-classical mapping of Sec. 12.31 Technically, the 
physics is contained in a model of the sine-Gordon type, where the action of a free field 
describing the phase dynamics is supplemented by a term arising from vortex tunneling 
events |3j- 

The model can be analyzed by conventional RG methods; interestingly, formally 
the same RG equations occur in the context of the Kondo model and will be discussed 
in Sec. 14.11 The RG analysis shows that the Kosterlitz-Thouless transition does 
not correspond to a critical, i.e., unstable fixed point in the RG sense, therefore 
"conventional" quantum critical behavior (as in Fig. and the corresponding power 
laws do not occur near in the vicinity of the transition point. Connected to that, 
the Kosterlitz-Thouless transition does not yield a singularity in any derivative of the 
thermodynamic potential at the transition, therefore it is sometimes called infinite-order 
transition. 

3.6.2. Fractionalization transitions. In the context of strongly correlated electron 
systems in dimensions d > 1, phases with "fractionalized" elementary excitations, 
carrying quantum numbers different from multiples of those of the electron, have 
attracted considerable attention in recent years. Fractionalization is usually discussed 
by re-casting the original model into a theory of new (fractionalized) quasiparticles, 
e.g., spinons and chargons, interacting with each other through a gauge field: the 
quasiparticles will carry a gauge charge. A fractionalized phase corresponds to a 
deconfined phase of the gauge theory, i.e., the effect of the gauge field is weak. In 
contrast, if the gauge field is confining, the fractionalized quasiparticles are bound 
at low energies, resulting in excitations with conventional quantum numbers. The 
fractionalization transition between such two phases corresponds to a confinement- 
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deconfinement transition of a certain gauge theory. The transition can be either driven 
by the condensation of topological defects of the gauge field, or by the condensation of 
matter (Higgs) fields which are coupled to the gauge field. 

The most prominent examples of gauge theories in the context of fractionalization 
are the ones with Z 2 [121 EH] and U(l) [Hj symmetry of the gauge field. Based on early 
work on the phase diagram of the U(l) gauge theory by Fradkin and Shenker [75], it is 
believed that Z 2 fractionalization can occur in d = 2, 3 whereas U(l) fractionalization 
can only occur in d = 3. The deconfined phases of those theories possess topological 
order which is associated with the suppression of topological defects in the gauge field; 
this order manifests itself also in a topological ground state degeneracy if the system is 
subjected to periodic boundary conditions. 

Fractionalization has been intensively discussed in the context of low-dimensional 
undoped and doped Mott insulators, with particular focus on high-T c superconductivity. 
In fact, Anderson's resonating valence bond state for two-dimensional antiferromagnets 
[7fi] is one of the first proposals for a fractionalized spin liquid. A powerful description 
of such a singlet liquid state is provided by Z 2 gauge theory models [73 [73] , see Ref . [2S] 
for further discussion. Experimental tests for fractionalization, based on the existence of 
vortices (so-called visons) in the Z 2 gauge field, have been undertaken, but no topological 
order has been detected so far in the cuprates. 

Significant progress has been made in constructing explicit models which show 
fractionalization; these involve simple models of bosons and fermions on regular lattices 
|74j . A notable achievement is the demonstration of fractionalization in the triangular 
lattice quantum dimer model [77j; furthermore indications for fractionalization have 
been found in numerical investigations of ring-exchange models on the triangular lattice, 
and charge fractionalization has been demonstrated in a correlated model on the 
pyrochlore lattice [75] . 

Recently, proposals for fractionalized phases in Kondo-lattice models have been put 
forward [3H]- Here, in the presence of strong quantum fluctuations the local moments 
can form a fractionalized spin liquid in a regime where the inter-moment interaction 
dominates over Kondo screening. The topological order protects this spin liquid against 
small perturbations, thus the conduction electrons will be essentially decoupled from the 
local moments. This peculiar paramagnetic phase has been termed "fractionalized Fermi 
liquid" , as it displays both Fermi-liquid-like particle-hole excitations and fractionalized 
spinon excitations, arising from the conduction electron and local moment subsystems, 
respectively. Increasing the Kondo coupling eventually drives the system through a 
(Higgs) confinement transition into the usual heavy Fermi liquid phase. Approaching 
from the Fermi liquid side, this zero-temperature transition can be understood as the 
breakdown of Kondo screening. The quantum critical region of this fractionalization 
transition displays novel critical behavior [SB] which may be related to experimental 
findings in heavy-fermion systems like CeCug-ajAu^ [3*5] . 

It is clear that fractionalization in d > 1 is a fascinating field which can lead 
to a plethora of new phenomena in condensed matter physics. Although experimental 
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evidence for the occurrence of such phases is lacking to date, we have to keep in mind that 
the experimental distinction can be very subtle due to the absence of a conventional order 
parameter; it could very well be the case that fractionalization is realized in recently 
studied strongly correlated materials. Therefore, further theoretical work predicting 
properties and proposing stringent tests of fractionalization is called for. 

3.7. Competing orders 

The physics of competing orders has been a central theme in low- dimensional correlated 
electron systems over the last years. In particular, in two dimensions extensive studies 
of microscopic models have shown that a variety of phases compete, and the actual 
ground state may depend sensitively on details of the model parameters. The interplay 
of different orders appears crucial for the understanding of the complex phase diagrams 
of various materials, with high-temperature superconductors being a paradigmatic 
example. 

On the theoretical side, competing orders imply the existence of multiple order 
parameters. A simple phenomenological description involves a Landau-type theory of 
coupled order parameter fields, in which the form of the allowed terms follows from 
symmetry arguments. Let us discuss the case of two order parameters <pA arid 4>b, which 
break different symmetries of the Hamiltonian: then the simplest allowed coupling is 
usually of the density-density type 1 1 2 1 0^ | 2 . The Landau theory will generically permit 
four phases: disordered, pure A, pure B, and coexistence of A and B. Tuning the system 
between the pure A and pure B phases leads either to two continuous transitions, with a 
coexistence phase or disordered phase in between, or to a first-order transition between 
A and B; the occurence of a single continuous transition requires additional fine tuning 
of parameters in the action. 

In the context of high-T c superconductors, the competition between magnetism and 
superconductivity has been studied extensively. On the one hand, static magnetism 
tends to suppress superconductivity, on the other hand, magnetic fluctuations are 
believed to mediate pairing (perhaps supplemented by phonons). The intimate 
relation between antiferromagnetism and <i-wave superconductivity has been suggested 
to arise from an (approximate) higher symmetry of the underlying model, namely 
an SO (5) symmetry unifying the two-component superconducting and the three- 
component magnetic order parameter [7JI]- Other proposals for orders competing and/or 
coexisting with superconductivity involve charge-stripe order [321 [33], nematic order 
[HQ] , circulating currents [2II1EI], and more exotic orders related, e.g., to spin-charge 
separation [73]. 

Recent neutron scattering and scanning tunneling microscopy (STM) experiments 
indicate that important competitors of superconductivity are incommensurate spin and 
associated charge density waves (CDW). The spin density waves (SDW) correspond 
to a collinear spin order at wavevectors Q sx = (tt ± e, tt) and Q SJ/ = (tt, tt ± e) where e 
depends on doping. On general symmetry grounds, such spin order (4> sa ) is accompanied 
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Figure 3. Left: Grayscale plot of the local DOS as measured by STM at the 
surface of r^S^CaC^Og+a; in a magnetic field of 7 T, integrated over an interval 
of subgap energies (after |H3)- The dark regions can be identified with the location 
of vortices. Around each vortex a checkerboard modulation is clearly visible; with the 
lattice constant being about 3.9 A for Bi2Sr2CaCu208+a;, the modulation period is 
found to be four lattice spacings. Right: Theoretical computation of the local DOS 
within a pinning model for the dynamic spin density wave fluctuations - here the 
unit of length is the lattice constant. (The crystal axes are rotated by approximately 
45 degrees in the left plot.) 

by charge order (0 C ) at wavevectors Q c = 2Q S ; formally this follows from the existence 
of a cubic term (0^ Q 0* + c.c.) in the field theory, which is allowed only if Q c = 2Q S due to 
momentum conservation. (Note that CDW order can exist without simultaneous SDW 
order.) Incommensurate spin and charge orders, coexisting with superconductivity at 
low T, have been detected in Nd and Eu-doped La 2 _ 1 Sr :r Cu04 jHOj. A number of other 
experimental results [%T] |H2j can also be explained by proximity to a quantum critical 
point at which spin/charge order disappears: even in optimally doped compounds which 
show no long-range spin order a suitable description of the collective spin and charge 
excitations is given by a theory assuming the vicinity of an ordering transition at lower 
doping. Then, the additional order parameters are fluctuating [23 EB], and will show 
up in suitable dynamic measurements. Charge fluctuations can also become static due 
to impurity pinning, and thus be detectable, e.g., by STM. An applied magnetic field 
can enhance competing spin/charge fluctuations or even stabilize static SDW/CDW 
order by suppressing the competing superconductivity due to vortex superflow |5H 125] . 
Recent STM experiments on optimally doped Bi 2 Sr 2 CaCu 2 8+x . in a magnetic field show 
that a checkerboard modulation occurs in the local density of states near vortex cores 
[52"] , see Fig. [3] This supports the idea of competing orders, and a likely explanation 
is that dynamic collective charge fluctuations are enhanced in the regions of weakened 
superconductivity and pinned by local imperfections, like the vortex cores themselves 
[EE]. A detailed modelling of the energy dependence of the STM signal indicates that 
it is best described by a checkerboard modulation in microscopic bond rather than site 
variables [SB] . 
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Figure 4. Schematic T — phase diagram (after Refs. ^2 E3 ) f° r the high 
temperature superconductors. The vertical axis denotes a parameter which can destroy 
antiferromagnetic order in the undoped limit, like the ratio of the near-neighbor 
exchange interactions; the horizontal axis is the hole concentration, 5. The ground 
state is insulating at small doping due to Coulomb interactions (shaded). All other 
phases at finite 8 show superconductivity which can coexist with bond/charge order 
and possibly collinear magnetic order at lower doping. The two arrows a) and b) 
correspond to the possible hole doping evolution of different compounds. 

At the microscopic level, the most important aspect of cuprates is that 
superconductivity is obtained from doping a Mott insulator. The generic phase diagram 
of doped Mott insulators is not understood to its full extent, however, it is clear that 
various spin and charge ordered phases, superconductivity, and more exotic phases could 
be realized. One simple scenario, based on large- N calculations for an extended t — J 
model, has been put forward in Ref . [H] : starting from a paramagnetic Mott insulator 
with bond order, doping naturally leads to superconducting phases coexisting with 
bond order, i.e., stripe-like charge density modulations. The resulting phase diagram is 
reproduced in Fig. |3J The interplay between superconductivity and bond order has been 
studied in some detail in Refs. jUHHH]: The ordering wavevector shows a characteristic 
evolution with doping, and period-4 structures have been found stable over a significant 
doping and parameter range. Phase diagrams with related physical ingredients, but 
some significant differences, appear in Refs. jHSl IHOj- We note that a number of 
experimental results indicate spin freezing into a glassy state below optimal doping, 
as has been pointed out by Panagopoulos et al. jSE] - this may be expected in the 
presence of disorder in a region of the phase diagram where spin ordering occurs in the 
clean limit. 

Besides the high-T c compounds, numerous other strongly correlated materials 
display competing orders. Evidence for co-existing antiferromagnetism and 
superconductivity has also been found in heavy-fermion compounds, e.g., CePd 2 Si 2 , 
CeCu2Si2, and Celn 3 . However, in these materials the nature of the superconducting 
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state has not been established to date. Other examples are colossal magnetoresistance 
manganites, ruthenium oxides, and organic conductors of the TMTSF and BEDT-TTF 
type. 

4. Boundary quantum phase transitions 

This section is devoted to the special class of boundary phase transitions. Such 
transitions occur in systems which can be divided into a bulk part with space dimension 
d and another part with dimension d^ < d - this "boundary" part can be either a surface 
or interface, or even a single impurity (c4 = 0) embedded into the bulk. At a boundary 
transition only the boundary degrees of freedom undergo a non-analytic change. At a 
continuous boundary phase transition, there is a singular part in the free energy which 
scales as L db where L is the linear dimension of the system. 

Boundary transitions have been extensively discussed in the context of surface 
transitions in magnets |87j . where under certain conditions the surface spins can order 
at a temperature above the bulk ordering temperature. The fascinating aspects of 
quantum criticality can also be found in boundary transitions. Here, we shall describe 
a few zero-temperature transitions in quantum impurity systems, which have been of 
much interest in diverse fields such as unconventional superconductors, heavy fermions, 
quantum dot systems, and quantum computing. 

4-1. Kondo effect in metals and pseudogap Fermi systems 

A magnetic impurity spin S embedded in a bulk system of non-interacting conduction 
electrons (Fig. Et) can be described by the so-called Kondo Hamiltonian, 

ff = 5> fc 4 <y c fcff + JxS-s(O) , (34) 

kcr 

where s(0) = Ekkw cl_ a a aa 'Ck' a ' is the conduction electron spin at the impurity site 
r = 0. The Kondo effect in metals, occuring at any positive value of the Kondo coupling 
Jk, is by now a well studied phenomenon in many-body physics. In essence, at low 
temperatures the conduction electrons and the impurity spin form a collective many- 
body state with zero total spin. The scattering of the conduction electrons off the 
impurity gives rise to non-trivial transport behavior, namely a temperature minimum 
in the electric resistivity of a metal with dilute magnetic impurities. For small Jk > 0, 
the low-energy physics of an isolated impurity is completely determined by a single 
energy scale, fceTx, where Tk is the Kondo temperature. The impurity spin is fully 
screened by the conduction electrons in the low-temperature limit, T <C T K . T K depends 
exponentially on the density of states (DOS) at the Fermi level, po = p( 6 f), and the 
Kondo coupling J K [88]. 

For ferromagnetic Jk < 0, the physics is different, and no screening occurs. 
Instead, in the low-temperature limit the impurity spin is effectively decoupled from the 
conduction band, and contributes a residual entropy of Sq = ks In 2. The two behaviors 
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Figure 5. a) Visualization of Kondo screening in metals: an impurity spin is 
surrounded by conduction electrons with primarily opposite spin, b) Kondo screening 
in superconductors is suppressed because the low-energy electrons are locked into 
Cooper pairs. However, for Jk larger than the Cooper pair binding energy screening 
is again possible - this can be interpreted as Bogoliubov quasiparticles bound locally 
to the impurity spin. 



can be understood in a weak-coupling renormalization group approach. If we allow for 
spin-anisotropic Kondo interaction, and define dimensionless couplings j± = PqJkl an d 
j|| = PqJk\\, the perturbative RG flow for small couplings takes the form (EE] 



sketched in Fig.|UK- There is a line of fixed points at j± = 0, jn < 0, and initial couplings 
with | J±\ < — J|| flow towards this line, whereas for all other initial values the system 
flows to strong coupling. The line \j±\ = represents a line of transitions of the 
Kosterlitz-Thouless type. Such a transition also occurs in the classical d = 2 XY model 
(Sec. IH.fj.ljl : it turns out that the RG equations for the Kondo problem are exactly 
identical to the ones for the XY model. In particular, no critical fixed point exists, and 
power law critical behavior is absent. 

This well-established picture of the Kondo effect has to be revised for systems 
with vanishing conduction band DOS at the Fermi energy. This is the case in hard- 
gap systems where p(e) = for energies smaller than the gap energy, |e| < A 
jHU EH], and in so-called pseudogap systems [HU E21 EH1 Ell with a power-law DOS 
p(e) = po\e/D\ r (r > 0). The former situation is realized, e.g., in semiconductors and 
s-wave superconductors, whereas the latter arises in semimetals and in systems with 
long-range order where the order parameter has nodes at the Fermi surface, e.g., p- and 
d-w&ve superconductors (r = 2 and 1). 

In systems with a hard gap in the conduction electron DOS, screening is impossible 
at small Kondo couplings Jk even at lowest temperatures, which is easily understood 
by the absence of low-energy fermions (Fig. EJd). Interestingly, in the particle- hole 
symmetric case screening does not occur even for large Jk, whereas in the absence 
of particle-hole symmetry an increase of the Kondo coupling eventually leads to a 
first-order boundary transition into a phase with Kondo screening - this transition is 
characterized by a simple level crossing between the doublet and singlet impurity states 



The case of a power-law DOS is even more interesting: the pseudogap Kondo model 
shows a continuous boundary quantum phase transition at a critical Kondo coupling, J c 




(35) 
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Figure 6. Renormalization group flow for the anisotropic Kondo model. The thick 
lines denote the separatrices corresponding to the phase transitions between local- 
moment (LM) and strong-coupling (SC) behavior, a) Metallic case, i.e., for a finite 
bath density of states po at the Fermi level. Here the transition is of Kosterlitz-Thouless 
type, b) Pseudogap case p(e) = po\e/D\ r (r > 0), showing a continuous transition with 
a critical fixed point CR. 

jUl E21 ESI • For small DOS exponents r, the physics is captured by the RG equations 

P(j±) = ~ r J±+j±j\\ , = ~ r 3\\ +3± , (36) 

see Fig. \§jp. There is now a critical fixed point at jj_ = j* = r. Initial values located 
below the separatrix flow to the local-moment fixed point (LM) at j = 0, describing an 
unscreened impurity, whereas large j flow to strong coupling (SC). The resulting phase 
diagram is shown in Fig. EK- The behavior at large r as well as the strong-coupling 
physics is not described by this RG; from comprehensive numerical studies 
based on Wilson's numerical renormalization group (NRG) approach [22] it is known 
that the fixed-point structure changes for r > 1/2, furthermore particle-hole asymmetry 
is a relevant parameter in the strong-coupling regime for r > (while being marginally 
irrelevant at r = 0). 

For small r, critical properties of the pseudogap Kondo model can be calculated 
using the weak-coupling RG sketched above. This can be understood as an expansion 
about the lower-critical "dimension" r = 0, similar to the expansion in (2+e) dimensions 
for the non-linear sigma model (|18|1. Very recently, it has become clear that r = 1 
plays the role of the upper-critical "dimension" in the pseudogap Kondo problem. The 
universal critical theory in the vicinity of r = 1 can be formulated as a crossing of 
singlet and doublet impurity levels, coupled to low-energy conduction electrons - this 
is equivalent to an infinite-?/ single-impurity Anderson model with pseudogap DOS 
[01]. For r < 1 the coupling between impurity and conduction band is relevant 
under RG transformations, and the model allows for a RG approach together with 
an expansion in (1— r), in analogy to the expansion in (4— e) dimensions of the 4 model 
(|17j) . Critical exponents can be calculated in a (1 — r) expansion, and hyperscaling is 
obeyed. Conversely, for r > 1 perturbation theory is sufficient, and the transition can 
be characterized as level crossing with perturbative corrections, where hyperscaling is 
violated [03]. We note that slave-boson mean- field theory [03] reproduces a transition, 
but does not yield sensible critical behavior. Alternatively, a dynamic multichannel 
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large-iV approach jHH] can be applied to the pseudogap Kondo problem, allowing for 
a complete analytic low-energy solution both in the stable phases and at the quantum 
critical points, with non-trivial critical behavior [HHj- 

Experimentally, signatures of non-trivial Kondo physics have been observed in 
impurity-doped cuprate superconductors. Nominally non-magnetic Zn impurities 
(replacing Cu) have been shown to induce quasi-free magnetic moments in their vicinity. 
The NMR data of Ref. [HZ] allowed to fit the local impurity susceptibility to a Curie- 
Weiss law, where the Weiss temperature (—0) can be roughly identified with the 
Kondo temperature. This experiment indicates a strongly doping- dependent Kondo 
temperature, with Tk in the superconducting state ranging from 40 K around optimal 
doping to practically zero for strongly underdoped samples, see Fig. ITJd. This qualitative 
change in the Kondo screening properties is a good candidate for a realization of 
the boundary quantum phase transition in the pseudogap Kondo model - there is 
no fundamental reason for this transition to coincide with any of the possible bulk 
phase transitions in the cuprates. As demonstrated in microscopic calculations 
the boundary transition from screened to unscreened dilute impurity moments in 
superconducting cuprates is likely driven by the increase of the superconducting gap 
upon underdoping and/or by increasing antiferromagnetic host spin fluctuations which 
suppress fermionic Kondo screening (Sec. 14.3)1 . Kondo screening of the Zn moments 
around optimal doping also leads to a characteristic low-energy peak in the tunneling 
signal as measured by STM, for further discussions see Refs. |§H1 HUUj . 

Kondo physics in a superconducting environment can be observed in quantum dot 
systems, too. A recent experiment |101j using a carbon nanotube dot coupled to Nb leads 
shows a sharp crossover in the transport properties as a function of T#/A, where Tk is 
the normal-state Kondo temperature of the dot and A the gap of the superconductor, 
consistent with the expectations within a hard-gap Kondo model. 

4-2. Spin-boson and Bose-Kondo models 

Systems of quantum impurities coupled to bosonic baths are an equally interesting model 
class; they have been first introduced in the context of the description of dissipative 
dynamics in quantum systems |l()2j . The simplest realization is the so-called spin- 
boson model, describing a spin or two-level system coupled to a single bath of harmonic 
oscillators. The Hamiltonian 



is a straightforward generalization of the spin-boson model to multiple baths: the a a are 
vector bosons and can be interpreted as spin-1 excitations of a "magnetic" bath. The 
couplings between spin and bosonic baths are completely specified by the bath spectral 
functions 




(37) 
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Figure 7. Left: Phase diagram of the pseudogap Kondo model. At T — the 
moment is screened at large Jk (SC), but unscreened at small Jk (LM). The two 
phases are separated by a continuous boundary transition at Jk — Jc] the quantum 
critical region is characterized by power laws, e.g., in the local susceptibility and 
the conduction electron T matrix 98, 95]. The crossover line between the quantum 
critical and the screened regime can be loosely associated to the Kondo temperature. 
Right: Temperature dependence of the inverse local susceptibility of dilute Li ions 
in YBa2Cu306+iE (after Ref. [57]). The circles (triangles) are data from optimally 
doped (underdoped) samples; the lines are Curie- Weiss fits. At optimal doping, the 
Weiss temperature changes significantly upon passing through the superconducting T Cl 
indicating that Kondo screening above T c is stronger than below T c . No change occurs 
in the underdoped sample, which should be attributed to the pseudogap above T c . The 
data below T c show a Weiss temperature (—0) of 40 K whereas (—6) is close to zero 
in the underdoped case - this suppression of Kondo screening likely corresponds to a 
boundary quantum phase transition. 

Of particular interest are power-law spectra 

J a {uj) = 2-n~i a u l ( T s u s , < oj < lu c , s > — 1 (39) 

where uj c is a cutoff, and the dimensionless parameters j a characterize the coupling or 
dissipation strength. 

In the conventional spin-boson model only j z ^ {pf x = j y = 0); then Eq. (|3*7|) 
describes a spin, tunneling between | f) and | j) via A, and being damped by the 
coupling to the oscillator bath. The particular case of s — 1 corresponds to the well- 
studied ohmic spin-boson model |102j . which shows a Kosterlitz-Thouless quantum 
transition, separating a localized phase at 7 > 7 C from a delocalized phase at 7 < 7 C 
102]. In the localized regime, the tunnel splitting between the two levels renormalizes 
to zero, i.e., the system gets trapped in one of the states | f) or | |), whereas the tunnel 
splitting stays finite in the delocalized phase. In the limit A <C uj c the transition occurs 
at 7 C = 1. 

Re-newed interest in spin-boson models arises in the field of quantum computation, 
for modelling the coupling of qubits to a noisy environment; here also the case of a 
subohmic bath appears physically relevant. It has recently been established, using both 
perturbative RG for small (1— s) and NRG, that the subohmic spin-boson model shows a 
continuous quantum transition between a localized and a delocalized phase for all bath 
exponents < s < 1 |l()3j . For fixed s there exists a single critical RG fixed point with 
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s-dependent exponents; as function of s these fixed points form a line which terminates 
in the Kosterlitz-Thouless transition point at s — 1. Interestingly, this behavior is 
somewhat similar to the one of the particle-hole symmetric pseudogap Kondo model 
jH2] described in Sec. 14. 1[ however, for s < 1 and r > the transitions in the two models 
are in different universality classes. 

Models of the form (|37J) with a vector bath describe impurity moments embedded 
into quantum magnets, and are sometimes called Bose-Kondo models. The bosonic 
baths represent the host spin excitations: in a quantum paramagnet their spectrum 
is gapped, with the gap approaching zero at a zero-temperature magnetic ordering 
transition. (In dimensions d < 3 the bath Hamiltonian has to be supplemented by 
bosonic self- interactions, as those are strongly relevant, see Sec. 12.71 ) At a bulk quantum 
critical point of an antiferromagnet in d = 2 + s space dimensions the bath spectra are 
given by Eq. ()39j) . For 1 < d < 3 the impurity model has remarkable properties 104J : 
the coupling between spin and bath is relevant in the RG sense and flows to an infrared- 
stable intermediate-coupling fixed point. The impurity spin is then characterized by 
power-law auto-correlations and various universal properties, e.g., a Curie response of a 
fractional effective spin. For a finite impurity concentration, this universal interaction 
between impurity moments and host spin fluctuations leads to universal impurity- 
induced damping of spin fluctuations in cuprate superconductors [104J - note that at 
the relevant spin fluctuation energies (e.g. 40 meV in YBa 2 Cus07) possible Kondo 
screening of the moments can be safely neglected, as experimentally Tk is 40 K or lower 

inn. 

4-3. Kondo effect and spin fluctuations 

In particular in the context of strongly correlated electron systems, which often feature 
Fermi-liquid quasiparticles and strong spin fluctuations at the same time, the question 
of the interplay between fermionic and bosonic "Kondo" physics arises. This naturally 
leads to so-called Fermi-Bose Kondo models where an impurity spin is coupled to both 
a fermionic and a bosonic bath. The renormalization group analysis shows that the two 
bath couplings compete, i.e., fermionic Kondo screening is suppressed by strong host 
spin fluctuations [SZ1 EHj ■ 

Recently it has been proposed that this interplay actually plays a role in 
high-T c cuprates, where NMR experiments indicate that dilute impurity moments are 
screened at optimal doping, but T K is essentially suppressed to zero in underdoped 
compounds |HZ|, see Sec. 14.11 A Fermi-Bose Kondo model, taking into account 
both the pseudogap density of states of the Bogoliubov quasiparticles and the 
strong antiferromagnetic fluctuations, provides a natural explanation: spin fluctuations 
increase with underdoping, thus strongly reducing Tk due to the vicinity to the boundary 
transition which exists even in the absence of a bosonic bath [99J. 

Interestingly, similar quantum impurity models also appear in variants of dynamical 
mean-field theories for lattice systems j^Sj. Motivated by neutron scattering experiments 
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[3~o] on the heavy-fermion compound CeCu 6 _ x Au x (see also Sec. I3.2.3|) . which indicate 
momentum-independent critical dynamics at an antiferromagnetic ordering transition, 
a self-consistent version of the Fermi-Bose Kondo model has been proposed to describe 
such critical behavior within an extended dynamical mean-field approach (EDMFT) 
|37j . In this scenario, the critical point of the lattice model is mapped onto a particular 
critical point of the impurity model. 

4-4- Multi-channel and multi-impurity models 

Kondo screening is strongly modified if two or more fermionic screening channels 
compete. Nozieres and Blandin |lU5j proposed a two-channel generalization of the 
Kondo model, which shows overscreening associated with an intermediate-coupling 
fixed point and non-Fermi liquid behaviour in various thermodynamic and transport 
properties. In general, such behavior occurs for any number of channels K > 1 coupled 
to a spin |, and does not require fine-tuning of the Kondo coupling. However, it is 
unstable w.r.t. a channel asymmetry; thus a two-channel fixed point can be understood 
as critical point between two (equivalent) stable single-channel fixed points. 

Many of the low-energy properties of the two-channel and related Kondo models 
have been studied using conformal field theory techniques (OBJ 11061 . Interestingly, the 
multi-channel Kondo fixed point is perturbatively accessible in the limit of large channel 
number (K ^> 1) pB1llU5j . Experimental realizations have been discussed in the context 
of rare-earth compounds |l()7j : furthermore proposals based on quantum-dot devices 
have been put forward. 

Models of two or more impurities offer a new ingredient, namely the exchange 
interaction, /, between the different impurity spins; it can arise both from direct 
exchange and from the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction mediated 
by the conduction electrons. This inter-impurity interaction competes with Kondo 
screening of the individual impurities; in Kondo lattice models it can lead to a magnetic 
ordering transition. The simplest model of two spin-| impurities has been thoroughly 
studied: here, a ground state singlet can be realized either by individual Kondo screening 
(if / < Tk) or by formation of an inter-impurity singlet (if I > Tr). It has been 
shown that these two parameter regimes are continuously connected (without any 
phase transition) as I is varied in the generic situation without particle-hole symmetry. 
Notably, in the particle-hole symmetric case one finds a transition associated with an 
unstable non- Fermi liquid fixed point [TUBl EHI • 

Quantum phase transitions generically occur in impurity models showing phases 
with different ground state spin. For two spin-- impurities, this can be realized by 
coupling to a single conduction band channel only. In this Kosterlitz-Thouless- 
type transition between a singlet and a doublet state occurs, associated with a secondary 
exponentially small energy scale in the Kondo regime |llUj . The physics becomes even 
richer if multi-channel physics is combined with multi-impurity physics - here, a variety 
of fixed points including such with local non-Fermi liquid behavior can be realized. 
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Experimentally, quantum dots provide an ideal laboratory to study systems of two 
or more "impurities" - note that the local impurity states can arise either from spin 
or from charge degrees of freedom on each quantum dot. In particular, a number of 
experiments have been performed on coupled quantum dot systems which can be directly 
mapped onto models of two Kondo or Anderson impurities, and some indications for 
phase transitions have been found In addition, experimental realizations of two- 

impurity models using magnetic adatoms on metallic surfaces appear possible. 

5. Conclusions and outlook 

This review has discussed aspects of zero temperature phase transitions in quantum 
systems. The introductory sections have highlighted the general properties of these 
quantum transitions, pointing out similarities and differences between classical thermal 
and quantum phase transitions. Subsequently, we have described a number of 
specific bulk quantum phase transitions in condensed matter systems, emphasizing the 
important role of low-energy fermionic excitations in determining the critical behavior. 
In particular, we have touched upon spin and charge ordering transitions, secondary 
pairing transitions in superconductors, as well as metal-insulator and superconductor- 
insulator transitions, and made contact to recent experiments in correlated electronic 
systems. Finally, we have discussed quantum transitions in impurity systems. Those 
transitions are a particular realization of boundary quantum phase transitions, where 
only a part of the system becomes critical. 

Let us recapitulate a few important aspects: The theoretical description of a 
particular phase transition usually starts out with the identification of the relevant 
variables, the most important one being the order parameter. Analytical studies then 
often proceed following the Landau- Ginzburg- Wilson approach: all degrees of freedom 
other than the order parameter fluctuations are integrated out, resulting in an effective 
theory in terms of the order parameter only; such a theory can often also be guessed 
using general symmetry arguments. However, in itinerant electron systems the existence 
of gapless particle-hole continua potentially leads to singularities in the order parameter 
theory, rendering it non-local. A more promising but technically difficult approach 
consists of not integrating out all degrees of freedom other than the order parameter, 
i.e., treating all soft modes in the system on the same footing. This has been followed 
for the ferromagnetic transition of both clean and dirty Fermi liquids [2I|, an d in the 
somewhat simpler situation of secondary pairing transitions in (i-wave superconductors 
|4*T] . A second caveat of analytical approaches is that they have to employ perturbative 
methods in analyzing the field theories; those methods can fail in strongly coupled or 
strongly disordered systems. Alternatively, high-accuracy numerical simulations either 
of simple microscopic models or of order parameter field theories can be used to access 
quantum critical behavior. 

Conceptually, quantum phase transitions open a field of fascinating physics, being 
connected to the peculiar properties of the quantum critical ground state. The absence 
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of quasiparticle excitations can lead to a variety of unsual finite-temperature phenomena 
in the quantum critical region, such as unconventional power laws and non-Fermi liquid 
behavior, which mask the properties of the stable phases. Recent years have seen an 
outburst of experimental activities studying quantum criticality in systems as diverse 
as high-temperature superconductors, quantum Hall systems, heavy-fermion materials, 
quantum magnets, and atomic gases. For explaining the plethora of experimental 
findings progress on the theory side is required, in particular in understanding dynamical 
and transport properties near quantum criticality, together with the influence of 
disorder. 

Quantum criticality can provide new perspectives in the study of correlated systems, 
where intermediate-coupling phenomena are hardly accessible by standard weak- or 
strong-coupling perturbative approaches. A promising route starts by identifying 
quantum critical points between stable phases, and then uses these as vantage points 
for exploring the whole phase diagram by expanding in the deviation from criticality. It 
is clear that we have just scratched the surface of much exciting progress to come. 
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